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ELECTRON  STATES  NEAR  SURFACES  OF  SOLIDS 


Daniel  C.  Mattis 

BELFER  GRADUATE  SCHOOL  OF  SCIENCE 
Yeshtva  University 
New  York,  New  York  10033 


ABSTRACT 


The  inatheniatical  theory  of  electron  eigenstates  near  the  surfaces 
of  solids  is  developed  in  stages.  We  first  study  the  eigenstates  of 
solids  which  terminate  at  a surface  but  are  otherwise  unperturbed,  and 
prove  that  near  the  surface  the  3-dimens ional  band  edges  are  "softened" 
and  van  Hove's  singularities  in  the  density  of  states  are  eliminated. 

A set  of  'Vacuum  states"  lying  primarily  outside  the  solid  is  constructed 
out  of  plane  waves  orthogonalized  to  the  eigenstates  of  the  solid.  This 
set  of  vacuum  states  is  not  orthononnal,  and  must  be  orthonormalized  by 
an  original  procedure,  different  from  that  of  Schmidt  (which  is  unsuit- 
ed). Effects  of  surface  perturbations  are  studied.  An  exact -me thoc 
is  elaborated  for  obtaining  eigenfunctions  in  closed  form,  for  a variety 
of  perturbing  potentials,  extending  arbitrary  distances  from  the  sur- 
face and  including  inter-band  matrix  elements.  It  consists  of  calcula- 
ting the  effects  of  one  surface  layer  at  a time,  and  cumulating  the  re- 
sults. The  intrinsic  instability  of  certain  surfaces  against  the  forma- 
tion of  bands  of  surface  states  is  shown  to  be  the  consequences  of  the 

vanishing  of  a bulk  quantity  a > one  of  the  components  of  the  inverse- 

22 

effective  mass  tensor,  at  certain  points  in  Che  B.Z. 


Resesrch  supported  in  part  bjr  the  Office  of  Naval  Research,  and  the 
National  Science  Foundation. 
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I . INTRODUCTION 

Vich  a vast  literature  now  appeared  on  various  aspects  of  the 
theory  of  surfaces  of  solids,  it  might  seem  that  the  scope  for  further 
investigation  is  thoroughly  circumscribed.  Tet,  even  a cursory  glance 
at  this  published  literature  reveals  the  lack  of  strong  mathematical 
underpinnings,  of  the  type  that  Floquet's  theory  and  group  theory  gave 
to  the  electron  states ( the  Bloch  functions^  of  the  periodic  solid. 

IIusSe  of  the  published  works  consist  principally  of  detailed  numerical 
investigations  of  the  Schrodlnger  differential  equation  in  the  neighbor- 
hood of  the  surfaces.  There,  perturbing  potentials,  both  one-body  and 
two-body,  are  taken  Into  account,  self-consistently  in  the  best  of 
these  investigations.  The  study  of  easily  solved  simple  models , such 
as  nearest -neighbor  simple  cubic  ("cubium")  structures,  yields  only 
qualitative  "seat  of  the  pants"  experience  but  no  generally  applicable 
formulas  or  theorems . 

In  the  present  work  we  have  taken  an  approach  which  we  hope  the 

Th- 
reader will  agree  to  be  justified  and  rewarding.  •Gar  idea  is  that  In 

any  surface  problem,  the  proper  basis  set  in  which  to  do  the  calcula- 
tions must  be  adapted  to  the  syrametry  of  the  problem,  in  this  case  the 

termination  of  the  solid.  If  this  feature  is  "built  in"  many  fea- 

\ 

tures  common  to  all  surfaces  can  be  Identified  in  advmce,  and  the  cal- 
culations for  specific  problems  reduced  to  a minimum. 

Then,  in  studying  the  effects  of  surface  potentials  on  theM  sur- 
face-adapted  functions,  we  find  ajs  now  quite  muted  and  pre- 

dictable; when  the  perturbing  potential  is  sufficiently  strong  or  the 
energy  contours  sufficiently  flat  (the  ■ 0 criterion^—  see  below). 
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there  eppeer  bound  states,  that  are  identified  as  surface  states.  In 
general,  the  basis  functions  are  scattered  among  themselves  by  the  per- 
turbations at  the  surface.  Calculations  of  the  scattering  and  of  the 
bound  state  properties  can  be  reduced  to  straightforward  simple  algebra, 
and  the  Introduction  of  more  sophisticated  effects  (many-body  corre- 
lations, principally)  then  Involve  no  more  difficult  concepts  or  calcu- 
lations than  for  the  bulk.  While  not  a universal  theory  of  surface  of 
solids,  the  present  work  does  outline  procedures  that  re-cycle  to  a 
maximum  extent  the  considerable  information,  already  available  about 
periodic  solids,  to  the  study  of  the  surface  related  properties.  Let 
US  now  examine  some  of  the  concepts  that  help  b*>i>g  jihn>*C . 


Although  chemists  describe  molecules  as  a collection  of  atoms  in- 
terconnected by  "chemical  bonds",  these  being  quantities  which  retain 
their  integrity  (length,  strength,  and  relative  angles)  in  a variety 
of  compounds,  little  use  of  this  important  concept  has  been 
made  in  solid-state  theory  outside  the  framework  of  such  approximate 
schemes  such  as  LCAO.  In  the  present  paper,  we  make  the  obvious  iden- 
tification of  the  bond  K(R^j)  connecting  an  atom  at  to  a second  atom 

at  R,  as  the  lattice  Fourier  transform  of  the  band  energy  c (k) , The 
J n 

exact  Bloch  energies  €j^(k)  and  corresponding  eigenstates  of  the 

petiovll.,  Icttlcc  cr?  Ingredients  in  the  present  formula- 


tion of  the  electronic  theory  of  surfaces  of  solids , 

(A)  Our  theory  of  surfaces  is  developed  in  several  stages . As  the 
first  and  most  Important  step,  bonds  which  connect  the  "inside"  atoms 
to  "outside"  atoms  are  cut  at  a surface  introduced  into  an  otherwise 
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•eml-tnf inlte  solid;  internal  bonds  as  well  as  the  internal  Haailltonian 
retain  their  bulk  values.  The  exact  new  eigenstates  are  then  alnply 
obtained  as  linear  combinations  of  the  Sloch  functions,  the  coefficients 
being  determined  by  a set  of  boundary  conditions.  While  this  procedure 
can  be  carried  out  l no 1 only  for  extremely  simple  examples  of 
nearest  neighbor  or  next-nearest  neighbor  bonding,  or  for  an  empty  lat- 
tice, our  method  is  so  formulated  that  in  general,  exact  eigenstates 
at  this  stage  can  be  easily  computed  using  simple  matrix 

algebra  on  the  known  bulk  solutions.  Moreover,  we  derive  a variety  of 
properties  of  the  eigenstates  of  the  terminated  solid  without  recourse 
to  numerical  methods.  For  example,  we  prove  in  all  generality: 

-There  are  no  "surface  states",  i.e.  no  bound  states  lying  in  the 
ranges  of  energy  forbidden  to  bulk  states.  The  energy  eigenvalues  of 
the  eigenfunctions  of  the  terminated  but  otherwise  unperturbed  solid 
interlace  the  Bloch  energies,  from  which  they  differ  only  to  order 

0(1/N^),  where  ■’number  of  planes  of  atoms  parallel  to  the  surface 

ee , 

-At  energies  c (k)  such  that  the  velocity  of  an  electron  perpendl- 
n 

cular  to  the  surface  v * A ^ 5e  Ck) /3k  vanishes,  so  does  the  relative 

z n z 

2 

probability  finding  it  at  a finite  distance  z from  the 

surface.  Thus  the  local  density  of  states  (l.d.o.s.)  at  planes  near 

ce  pifiCk  are  lacking  the  accumulated  contributions  from  band 

edges  and  saddle  points  in  the  energy  spectrum,  with  consequent  erosion 

of  the  bend  edges  and  elimination  of  the  van  Hove  singularities  charac- 

3 

terlstic  of  the  bulk  d.o.s.  and  band  structure. 
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(B)  Turning  to  the  second  stege  of  the  theory,  the  instabilities  of 


the  new  eigenfunctions  against  simple  model  surface  perturbations  is 

examined.  Generally,  surface  potentials,  or  changes  in  the  bond 

strengths,  can  "pull"  eigenstates  into  the  surface  or  "push"  them  out 

of  it.  More  importantly,  such  perturbations  can  cause  bound  states, 

i.e.  surface  states,  to  appear  -oesaidt^the  continuum  of  the  Bloch 

states.  Thus,  for  a given  point  k ,k  in  the  two-dimensional  (2D) 

X y 

Brillouin  Zone  (B.Z.)  of  the  surface,  we  could  have  states  lying  out- 
side the  continuum  of  energies  e (k  ,k  :k  ) (defined  by  allowing  k a to 
wander  over  the  Interval  -r:  to  provided  the  perturbation  is  suffi- 
ciently strong  or  the  incipient  instability  sufficiently  great.  At 


this  point  we  show  that  the  instabilities  arise  at  or  near  certain 


unusual  points  in  the  2D  B.Z.  where  the  inverse-effective -mass  parameter 
-2  2 2 

o^^(k^,ky)  s A 3 c^(k)/Bk^  vanishes,  the  derivative  being  evaluated 

at  the  values  of  k for  which  e (k  ,k  ;k  ) has  its  maximal  and/or  mini- 
2 n X y 2 

mal  values,  at  fixed  n,k^  and  k^.  Only  a knowledge  of  the  bulk  band 
structure  and  a choice  of  surface  orientation  relative  to  the  crystal 
axes  (to  determine  the  orientation  of  the  2-axis,  normal  to  the  sur- 
face) are  required  to  locate  the  xeros  of  *nd  thus  to  establish 

the  bands  and  surface  orientations  most  likely  to  produce  surface  states ^ 
(C)  In  the  third  stage,  more  general  perturbations  are  considered 
and  the  mathematical  procedure  for  obtaining  the  new  eigenstates  and 
eigenvalues  — both  in  the  continuum  and  for  the  bound  surface  states 
— is  developed.  It  is  only  for  large  perturbations  that  we  can  expect 
to  see  surface  states  lying  within  the  bulk  band-gaps,  as  weak  per- 
turbations produce  surface  states  whose  energies  usually  overlap  parts 


4 


of  the  bulk  continuum. 


■do  not  mix  as  long  as  they  belong  to 


different  values  of  n or  k and  k .f  The  various  contributions  to  the 

X V 

perturbing  Hamiltonian  Include:  a changed  Madelung  potential  (zero  out- 
side the  solid,  varying  progressively  to  its  bulk  limiting  value  at 
large  z) , a macroscopic  dipole  moment  due  to  the  imbalance  of  forward- 
backward  charge  on  orbitals  of  surface  atoms,  relating  to  changes  to 
the  usual  Hartree-Fock  and  correlation  energies  due  to  termination  of 

4 

the  crystal.  In  addition  to  scattering  of  the  various  eigenstates, 

such  perturbations  can  also  mix  different  bands  (conserv’ing  k ,k  how- 

X y 

ever),  resulting  in  a distortion  (rotation,  elongation,  or  whatever) 
of  the  original  bonds.  They  can  also  admix  the  "vacuum  states",  ori- 
ginally decoupled  by  the  breaking  of  surface  bonds. 

For  the  purposes  of  a realistic  calculation,  one  must  therefore  re- 
introduce a complete  set  of  states  for  the  outside  space  — the  "vacuum 
levels".  These  may  be  chosen  as  a complete  set  of  plane  wave  states 
displaced  by  an  energy  determined  by  the  magnitude  of  the  work  function, 
and  orthogonalized  to  the  eigenstates  of  the  solid  near  the  surface. 

A procedure  for  doing  this  is  outlined.  The  perturbations  described 
in  the  previous  paragraph  have  matrix  elements  connecting  the  inside 

states  to  the  vacuum  states  as  well  as  to  each  other,  and  so  we  had  to 

5 

have  a scheme  far  calculating  them. 

This  work  is  intended  to  provide  a thorough  Introduction  to  a va- 
riety of  new  mathematical  aspects  only;  we  Intend  to  present  applica- 
tions in  subsequent  publications.  Problems  concerning  "reconstruction" 
of  surfaces  are  also  deferred. 
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II,  SOME  FAMILIAR  BIT  USEFUL  CONCEPTS 

II 

The  bulk  Schrodtnger  equation  Hi (r)  » ci (r)  decenntnes  the  complete, 
orthonormal  set  of  Bloch  states  ♦ . These  can  be  written  in  two  distinct 

ways : 


♦ ("r)  » N ‘ u (r)l  or 

nk  x,y,2  — nk 


(2.1) 


, l(k-R  ) 

N ^ He  ^0  (r-Rj, 
x,y,z  n'  y 

i 


in  which  u^j^(r)  is  a function  having  periodicity  of  the  lattice,  and 

0 (r-R.)  is  the  Wannier  function  for  the  nth  band  at  the  R.th  lattice 
n'  3'  3 

site.  The  orthonormal  Wannier  function  is  localized  within  a small 
neighborhood  of  its  nominal  site,  decaying  exponentially  at  large  dis- 
tances from  it.  The  energy  eigenvalue  of  the  above  Bloch  state  is: 


' N 


1 


ik-R.  . 


Tie 


x,y,z  R^  Rj 


(2.2) 


It  is.  In  the  extended  zone  scheme,  periodically  extensible  into  the 
second,  third,  etc.  Brlllouin  Zones.  For  our  purposes,  knowledge  of 
(2.2)  for  k in  the  first  Brillouin  Zone  is  sufficient.  For  the  study 
of  surfaces,  an  even  smaller,  two-dimensional,  Brillouin  Zone  is  impor- 
tant. Let  the  normal  to  the  surface  plane  define  the  z-direction,  so 
that  the  surface  defines  the  x-y  plane.  The  new  boundary  conditions 
at  the  surface  and  the  perturbing  potentials  localized  near  the  sur- 
face may  or  may  not  affect  the  translational  symmetry  in  the  x-y  plane. 
When  they  do,  we  say  the  surface  is  "reconstructed".^  In  the  present 
work  we  avoid  these  complications.  We  shall  always  assume  the  aurface 
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to  be  unreconstructed,  l.e.  we  shall  make  strong  use  of  the  hypothesis 
that  the  translational  synmetry  In  the  x-y  plane  is  maintained  exactly 
«s  in  the  bulk  even  though  this  symmetry  is  destroyed  in  the  t-dlrec* 
tlon  by  the  surface. 

At  every  fixed  value  of  k^  e (k^,k^)  the  energy  e(kjj,k^)  achieves 

a minimum  and  a maximum,  at  k , and  k respectively.  These  will 

X min  z max 

be  denoted  the  "relative  band  edges"  to  distinguish  them  from  the  "ab- 
solute band  edges"  at  which  the  entire  nth  band  achieves  its  lowest  and 


min 


highest  energies.  We  denote  the  lower  relative  band  edge  (kjj)  a 


nd 


. max,,  s 

the  upper  one 


J^t  any  k in  the  nth  band  one  defines  an  inverse-effective-mass 
a 

tensor  with  nine  components  given  by: 


(2.3) 


In  surface  studies,  one  is  particularly  concerned  with  just  one  of 

these  components,  a , evaluated  at  k , and  k . We  Introduce 
’ n,rz  z min  z max 

the  self-explanatory  notation  (kn)  and  c?”**  (kn) , both  functions 

^ n,zz2  H n,x2'  Ir 

only  of  kjj;  quantities  which,  together  with  E°^”(kjj)  and  £*°*(^|j) 
must  evaluate  using  knowledge  of  the  usual  (bulk-)band  structure. 

To  fix  ideas,  consider  the  hypothetical  simple-cubic  tight-binding 
"cublum",  which  has  a band  structure: 


c “ K_(cos  k a -f  cos  k a + cos  k a ) 
n “ x o y o y o 


(2.4a) 
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for  the  surface  perpendicular  to  one  of  the  crystal  axes,  and 


i 

I 

1 


e “ K (cos  k a -+  2 cos  k a 2 ^ cos  k a 2 (2.4b) 

n n xo  yo  zo 


for  a surface  in  the  (0,1,1)  direction.  For  the  first,  (0,0,1)  orient- 
ation, ve  obtain  assuming  K < 0, 


k 


2 nin 


0, 


k 

z max 


, min, max 
TT/a,  e “ 


K (cos  k a -(-  cos  k a ) 4.  K 
n X o yon 


, min  ,max 
and  or 

zz 


± iKjj!  (setting  t 


1). 


whereas  for  the  second,  (0,1,1)  orientation. 


k , • 0 for  cos(  k a 2 *)  > 0 and  k “ Tr2^/a  for  cos  (k  a 2 0, 

z min  y o z min  o ' y o ' ’ 

i. 

with  k 


■=  t.2~/b  in  the  former  and  0 in  the  latter  case.  Then, 
z max  o ’ 


min  ,max 


e 

and 


K (cos  k a i 2 I cos  k a 2 ^1) 
n X o ' y o ' ' 


min, max  . l ,,  2 , 

a “ i K a cos  k a 2 ^ 

22  n y o ' 


(2.5) 


For  the  first  orientation,  the  two-dimensional  Brillouin  Zone  (2DBZ) 

is  a square  extending  from  -tt  to  -fTt  for  both  k a and  k a . For  the 

X o y o 

second,  it  is  a rectangle,  with  extending  from  -n  to +rt  while 

Vi?  tain 

extends  from  -tt2  to-¥‘n2‘.  We  display  contours  of  constant  e and 

the  behavior  of  in  Figures  1 and  2 (e  and  or”®*  are  similar). 

Z2  max 

In  equations  (2.1)  and  (2.2)  the  factors  N refer  to  the  total 

number  of  cells  in  the  3D  crystal.  With  the  total  number  of  planes 


i 


1 


parallel  to  the  aurface  plane,  and  N the  number  of  cells  In  each 

^ I y 

plane,  we  have  obviously: 


N - N N (2.6) 

x,y,z  z x,y  ' 

K^(R^^)  in  equation  (2.2)  is  the  "bond  .crrength"  or  Hamiltonian 

matrix -element  connecting  the  Wannier  function  in  the  nth  band  in  the 

cell  centered  at  R,  to  that  at  R..  While  in  the  periodic  solid  K is  a 
i j n 

function  only  of  the  relative  vector  R. . e R.-R  , introduction  of  a 

ij  1 J’ 

surface  complicates  matters.  Zero  Boundary  Conditior  ■ reouire  that  we 

cut  all  bonds  that  extend  outside  the  solid,  i.e.  rea-'ire  us  to  set  all 

K “ 0 for  which  either  R.  or  R.  lie  outside  the  solid,  jfhis  procedure 
n i J 

will  be  examined  in  the  following  Section. 

In  the  tight -bind  ing  scheme'^  the  bond-strengths  K (R , ) are  esti- 

n ij 

mated  first,  and  the  band -structure  calculated  therefrom.  In  more  mod- 
8 

em  schemes,  the  band  structure  (R)  is  obtained  by  direct  solution 

of  Schrodinger ' s differential  equation  appropriate  to  the  periodic 

solid,  with  or  without  correlation  and  exchange  corrections.  Once  this 

is  done,  we  can  convert  (2.2)  to  obtain  the  K 's : 

n 


-ik-R 

I 

l^t^B.Z. 


c^(k)  - 


In  the  Wannier-functlon  representation,  the  eigenstates  of  the  Hamilton- 
ian for  electrons  in  the  solid  take  the  general  form: 


*(r)  - E U(R,  )0^(r-R  ) 
0 J n J 
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where  the  V's  are  the  cotrponents  of  the  eigenvectors  of  the  N x N tns- 
trix  (2.7),  and  are  to  be  normalized  to  unit  length.  By  using  a trans- 
lation-operator representation  of  (2.7)  we  can  obtain  a Schrodinger- 
like  equation  for  these  amplitudes: 


H U(R)  = E U(R)  (2.9) 

n 

in  which  the  "Hamiltonian"  is  a generalized  translation  operator: 

H s € (-i-s/cR)  -IK  (RJ  exp(R.-S/aR)  (2.10) 

n n - n j j 

j 

When  subject  to  PBC  the  normalized  solutions  of  (2.9)  are  evidently 
plane  waves  N ^ exp(ik-R.)  and  (2.8)  reduces  to  the  Bloch  form  (2.1), 

^ >y  * 2 J 

with  the  energy  eigenvalue  E = e (k)  . 

n 

In  tight-binding,  the  Wannier  functions  are  presumed  known,  or 
are  approximated  by  the  corresponding  atomic  orbitals.^  In  the  more 
modem  procedures,  the  Bloch  states  or  their  periodic  components  u , (r) 

~ K 

are  obtained  by  a direct  computer  solution  of  the  appropriate  Schrodinger 

g 

differential  equation;  we  may  then  use  these  known  eigenfunctions  to 
construct  the  Wannier  functions,  by  inverting  (2.1): 


0 (r-RJ  = n"^ 

n J JC.y.2 


-ik-R 


k C l^'^B.Z. 


3 lk-(r-R  ) 


(2.11) 


We  shall  henceforth  assume  these  functions  to  be  kncm. 
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Ill,  ZERO  BOUNDARY  CONDITIONS 


A terminated  solid  consists  of  atomic  planes  at  Z.  » ni.a,  with 

J J 

TDj  •»  1,2,. ..,N^.  The  interplanar  distance  £ is  a function  of  the  ori- 
entation of  the  crystal  axes,  of  course  Ca  • a for  the  (001)  example 

o 

in  Section  II  and  a * a^2  ^ in  the  second,  (Oil),  example).  The 

existence  of  a surface  at  « 0 ellmins  tes  the  PBC  in  favor  of  re- 
quirements  that  the  amplitudes  U (defined  in  the  preceeding  Section) 
vanish  at  m^  * 0 , -1 , -2 , . - . . These  are  denoted  the  zero  boundary 
conditions  ZBC. 

We  take  advantage  of  periodicity  in  the  x-y  plane  to  factor  the 
amplitudes ; 


U . (R.)  = N ^ exp  i(kX.-(-  kY.)  U ,(Z.) 
nk  j'  x,y  ^ " X J y y nk''  j'' 


(3.1) 


where  the  U(Z) 's  satisfy  the  difference  equation: 


(3.?1 


for  which  the  translation  operator  is  given  as; 


V m 


(3.3) 


Alternately , 


V 


■+?  paS/Sz 

Z K(p)  e 

S--P 


(3.4) 


in  which 


11 


e 


(3.5) 


-fn/«  -ipak 

— n/a 


by  obvious  special Iratlons  of  the  preceeding  Section.  It  is  important 
to  emphasize  that  the  coefficients  K(p)  will  vary  with  planar  wave-vec- 
tor kit  E (k  ,k  ) and  band  index  n and  should  be  so  labelled;  K ,,,,(?)• 

11  ' x’  y'  - n.k'ii 

For  all  reasonable  band  structures,  one  may  without  significant  error 
cut  off  the  sum  in  (3.4)  beyond  a finite  maximum  integer  P,  which  rare- 
ly exceeds  0(10)  but  may  often  be  as  small  as  1 or  2 . We  assume  the 
existence  of  such  a cut-off  P,  and  amplify  much  of  the  material  in 
refere.'ce  1 for  the  sake  of  a complete  presentation. 

We  start  with  the  inverse  of  (3.5): 


ik  a 

c(k  ) - K(p)  (e  ^ (3.6) 

z p=-P 

Here  and  elsewhere  in  this  paper,  we  suppress  both  the  band  index  n_  and 
the  planar  wave-vector  kj^  for  reasons  of  typographical  clarity.  VH:iere- 
ever  this  could  occasion  some  confusion,  we  restore  all  the  parameters. 
We  also  use  the  notation  t ■ (k||,n)  on  occasion,  e.g.  K^(p)  is  the 
quantity  defined  in  (3.5)  and  E^(k^)  in  (3.6),  the  latter  also  denoted 

c„(k). 

The  equation  (3.6)  will  now  be  viewed  as  a 2Pth  order  equation  for 
exp(ik^a),  having  in  general,  2P  distinct  complex  roots  . We  label  these 
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* expCik^''  vlth  index  o«  1,2,...,2P.  The  most  general  solu- 

tion of  (3.2)  at  a fixed  E “ e(k  ) for  a real  specified  value  of  k is 


U(2j) 


2P  (a)  2P 

Z A exp(ik  '■  ''m.a)  « I A C ^ 
Qts\  a z j Qcitt'a 


(3.7) 


(or)  , 


6.  e.g. 


in  which  the  given  k is  one  of  the  members  of  the  set  of  k 
2 2 

k “ k The  coefficients  A are  now  obtained,  under  the  simplify- 

z z ct  ^ 

ing  assumptions  of  no-spin-orbit  coupling,  which  allows  us  to  assume 
the  K's  in  (3.4)  - (3.6)  real  and  symmetric  under  the  interchange 

P — ”P  • 


(i)  By  an  obvious  synsnetry,  for  every  root  k there  is  a second 
root  -k.  We  know  one  real  root,  k = the  specified  k, . Including  its 

Z z 

■negative,  there  is  a total  of  2Q  real  roots,  with  Q > 1.  The  remain- 
ing 2(P-Q)  roots  are  complex,  so  half  of  them  (those  with  negative  im- 
aginary parts)  must  be  discarded  because  they  lead  to  growing  (and 
therefore  un -normalizable)  exponentials  inside  the  solid  (m^  > 0).  The 
magnitudes  are:  | c!  * 1 for  2Q  real-k^  roots  and  ! C!  < 1 for  the  re- 
maining P-<5  complex-k^  roots  that  are  retained. 

(ii)  An  alternate  formulation  of  (3.6)  regards  this  as  a Pth  order 
polynomial  equation  in  cos(k^a).  Because  of  the  real  coefficients,  the 
roots  come  in  complex  pairs  such  as: 
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(3.8) 


{cos  k a - u + lw)(cos  k a - u - Iw)  = 0 
z z 

If  k is  here  an  acceptable  root  (l.e.  In  k > 0),  then,  by  in- 

Z 2 

spection,  so  is  and  the  two  are  distinct  unless  they  are  purely 

laaginary.  Summarizing:  the 

(a)  purely  real  's  come  in  i pairs,  the 

(b)  complex  's  come  in  pairs:  k^^®^  and 

(c)  Pure  imaginary  are  necessarily  unique. 

(iii)  We  now  impose  2BC  on  (3.7),  setting  U(Zj)  * 0 for 

m^  * 0,-1, . . . ,-(P-l) . There  is  no  requirement  at  m^  < ~P  as  the  bonds 
within  the  solid,  which  are  to  be  ’'broken"  by  the  imposition  of  the 
boundary  conditions,  extend  out  only  as  far  as  -(P-1).  Thus  the  P + Q 
non-zero  coefficients  in  the  expansion  (3.7)  are  constrained  by  a 
set  of  P linear  homogeneous  equations,  and  are  reduced  to  a set  of  at 
most  Q indep>endent  members.  The  net  number  Q is  determinable  as 

follows:  we  study  the  behavior  of  e (k  ) as  k is  varied  in  the  real 

z z 

interval  0 < k^  < rt/a.  In  the  range  where  e is  single-valued,  Q “ 1. 

In  the  range  (if  any)  where  c is  double-valued,  Q •=  2 , etc.  This  is 
illustrated  in  Figure  3.  We  note  that  the  net  number  of  solutions  is 
precisely  half  of  that  obtained  using  PBC,  a not  unexpected  result. 

(iv)  Orthogonalization  and  Normalization  of  solutions.  We  dis- 
tinguish two  types  of  sums , those  which  depend  on  N : 

z 

N 

IC  N . only  for  k real.  i.e.  IC  I - 1 (3.9) 

O'  z ~ ^ z ' *or' 

and  those  which,  in  the  thermodynamic  limit  *.  independent  of  N^ : 
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(3.10) 


N,  * 

««  'b’ 


08  IT  _ T 


'-C^e 


The  latter  occur  when  either  or  complex  or  inaginary, 

or,  in  general,  for  a ¥ ^.  (We  recall  that  the  imaginary  parts  are  al- 

ways  positive;  for  real  we  may  here  (but  not  in  (3.9))  introduce 

an  infinitesimal  positive  imaginary  pert  to  speed  the  convergence.) 

(a)  Let  us  now  normalire  a function  of  the  type  (3.7)  when  Q = 1. 

There  are  P+1  terms  in  the  expansion  (3.7):  the  specified  k and  -k  , 

z z 

a nertain  number  S of  purely  imaginary  and  the  remainder 

(P-l'S)  complex  pairs  of  the  type  and  Thus,  (P-l-S)  is 

an  even  integer  or  zero.  We  write  U as  follows: 

U(Zj)  - 

The  ratios  a^  * A^/Aj^,  a^  * A^/A^,...  are  determined  by  the  P boundary 

conditions.  We  have  chosen  k « k , k » — k , so  that 

z z’  z z 

j^l)  - lC2i“  1 other  ) ^1  < 1.  In  part  (b)  below,  we  shall 


prove  Aj  ¥ 0. 

Because  both  the  Hamiltonian  and  the  boundary  conditions  are  real 
the  solutions  (3.7)  are  real  (or  can  be  made  real  by  a trivial  (constant) 
multiplicative  factor).  Thus,  a2^l“  ^i*''  coefficients  a^A^  for  the  S 
terms  involving  the  pure  imaginary  k^^°”^'s  are  all  real;  and  the  coef- 
ficients a , a , of  the  (P-l-S)  conjugate  pairs  are  themselves  complex 
n nf  1 

conjugates:  . 

The  normalization  requirement  is  thus 
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(3.12) 


i 


1 

I 


- [2N^  + 0(1)1 

by  use  of  (3.9)  end  (3.10)  and  |a2]^  “ 1.  Thus,  is  a constant: 


Ia,|  - (1/2N  ) 
1 z 


(3.13) 


independent  of  e(k^)  and  of  a^^  for  i > 3,  a fortunate  s itr.pl  if  icat  ion . 

(b)  The  degenerate  cases  for  Q > 2 are  treated  analogously,  by 

constructing  Q linearly  independent  functions  of  the  type  (3.11)  each 

containing  one  or  more  of  the  real  k 's  (and  their  mates  -k  ),  ortho- 

z z 

gonalizing  them,  and  then  normalizing  each. 

This  requires  us  to  prove  that  one  cannot  construct  a solution 
utilizing  only  the  P complex  or  imaginary  allowed  Such  a proof 

is  also  required  in  part  (a)  preceeding,  to  justify  the  assumption 
¥ 0,  on  which  the  entire  analysis  (3.11)  - (3.13)  was  based. 

Supposingone  could  construct  a solution  U using  anv  subset  (P'  < P)aaP' 
distinct  k^'  '^'s,  labeled  a“  Qfl  ,Qf  2 , . . . ,Qf  P , such  as: 

/ 

«■“> 

with  the  ZBC  applied  at  m » 0 ,-l , . . . ,-(P-l)  to  determine  the  A^'s. 

These  P^equations  may  be  conveniently  auimnarized  in  matrix  form: 
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T — 


■ 

— — 

A, 

8 

1 

^2 

• 

cs 

• 

M 

• 

A . 

P' 

u - 

“ — 

in  which  the  elements  of  the  PxP  matrix  are: 

1-s  * 

• Cj.  » column  index,  r - 1,...,P 

and  row  index  s — 1,...,P^ 


(3.15) 


(3.16) 


Equation  (3.15)  can  be  satisfied  Iff  the  determinant  of  M vanishes.  This 
is  recognized  as  Vandermonde's  determinant,  having  the  value: 


Det  1Im!| 


i "V 


<3.17) 


q>r  T*1 

which  cannot  vanish  as  the  Cj.'s  are  all  distinct,  (J.E.D. 

The  above  has  as  one  important  consequence,  that  no  decaying  solu- 
tion (alias  surface-state)  can  be  found  to  the  unperturbed  ZBC  problem, 
i.e.  all  states  in  this  model  have  at  least  one  component  with  I d * 1 
that  survives  in  the  bulk.  A second  corollary  is  that  no  new  states 
are  created  in  the  energv-gap  region  by  the  Imposition  of  ZBC,  as  Q “ 0 
there . 

The  next  step  is  to  construct  Q distinct  functions  of  type  (3.11), 
each  leading  off  with  a different  chosen  from  the  set  of  Q real, 
positive,  specified  k^'s  illustrated  in  Figure  3.  The  reals  again  come 
in  pairs:  thusla^l  • 1;  however  we  can  no  longer  assume  that  ■tid 


I. 
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the  succeeding  ones  are  complex,  as  a number  of  them  may  belong  to  the 
set  of  reals.  Let  denote  the  set  of  reals  contained  In  one  of 

the  linearly  independent  solutions;  it  contains  at  least  2 members,  and 
at  most  2Q.  Let  be  the  set  of  reals  in  a second  solution.  The  sets 
and  R^  may  be  distinct,  in  which  case,  by  equations  (3.9)  and  (3.10) 
the  two  solutions  are  effectively  orthogonal  (in  the  limit  “)  • 

However,  they  may  have  one  or  more  pairs  of  reals  in  common,  as  in  the 
following  example : 


A(i 


m 


a. 


+ a. 


and 


02-  B(C^i  ...). 

the  dots  indicating  remaining  terms  all  involving  complex  or  imaginary 
k^'s.  The  normalization  requirement  for  L'2  leads  to; 

|b!^  - [2N  1'^ 

'■  z 

We  Schmidt-orthogonallze  the  two  solutions,  replacing  Uj^  by: 

- (N^A*B(a*  i «3))U2 


The  solution  can  now  also  be  normalized.  Note  that  the 
Schmidt-orthogonalization  and  the  normalizations  are  all  to  be  carried 
out  without  regard  to  evanescent  complex  or  Imaginary  k 's.  This 
may  be  regarded  as  one  of  the  important  simplifications  afforded  by  pro- 
ceeding to  the  infinite  "thermodynamic"  limit,  and  distinguishes  the 


3 


1 

. 

3 


i 
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present  procedure  from  those  In  which  the  solid  is  modeled  by  a finite 
■number  of  parallel  planes. 


(v)  Deep  in  the  bulk  (m^  — »)  the  evanescent  components  disappear, 


so  for  any  properly  normalized  ^nd.  If  need  be,  Schmidt-orthogonalized) 
solution  U (to) : 


lim„  lu(m)r  = 1/N 

m — * ‘ ' 2 


(3. IS) 


at  any  energy  c(k^),  for  any  degree  Q of  degeneracy.  We  define  the  various 
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local  density-of-ststes  functions  for  the  mth  plane  in  the  nth  band  as; 


,-l 


O'.® 


(3.19) 


By  the  preceeding,  we  see  that  the  bulk  limit  pj(E)  is  attained; 


.-1 


■m 


li®,  e P^(E)  - I 6(E-€^(k||,kp) 


(3.20) 


which  includes  all  degeneracies  in  the  nth  band.  The  summed  density-of- 


states  p(E)  e Lp^(E)  is  sometimes  also  of  interest.  The  evaluation  of 


such  quantities  at  or  near  the  surface,  where  the  evanescent  contribu- 
tions are  in  principle  important,  is  the  topic  of  Section  5. 


IV.  SURFACE  AMPLITUDES  (UNPERTURBED) 

3 7 11 

The  nearest-neighbor  "tightest-binding"  (P“l)  ■model  * is  readily 


11 


solved  by  the  above  procedures,  to  obtain 


I* 

i 
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(2/N^)‘  sln(k^ 


The  bulk  limits  (3.18)  - (3.20)  can  be  explicitly  verified,  and  the  be- 
havior near  the  surface  can  be  determined.  In  particular,  ve  aee  that 
at  finite  distance  m from  the  surface  the  amplitudes  near  the  band  edges 
(at  k 0 and  rr/a)  vanish  as  (k  - k . ) and  (k  - k ) respec- 
tive ly . 

In  this  Section  we  shall  prove,  confirming  our  earlier  conjecture,^ 
that  this  "erosion"  of  the  band  edges  is  a perfectly  general  phenomenon. 
We  shall  then  analyze  the  behavior  of  the  surface  density-of-states 

•it 

function  and  of  such  Important  quantities  as  l!j^(m)Uj^,  (m  ) near  the 
relative  band  edges,  as  required  for  the  analysis  of  the  perturbed 
surfaces  in  the  following  Sections. 


For  the  sake  of  definiteness,  we  study  only  the  band  minima;  the 
maxima  are  studied  in  precisely  the  same  way,  and  do  not  require  eudi- 
tlonal  commentary.  Note  in  passing  that  a relative  band  extremum  can 
be  non -degenerate  (Q=l)  only  if  it  occurs  at  k^  ••  0 or  -n/a  as  in  (4.1) 
above.  If  it  occurs  at  an  arbitrary  point  within  the  Zone,  Q > 2;  a 
local  minimum  such  as  point  ( kj^ ) in  Figure  3 has  Q > 3. 

(a)  To  start,  let  Q * 1,  P arbitrary.  We  use  equations  (3.11)  - 
(3.13)  for  the  form  of  Uj^(m) . We  now  write  the  P "ZBC  equations"  for 

m “ 0,-1, , plus  one  equation  at  m ■!  for  Uj^(l)  , extending  the  matrix 

formulation  introduced  in  equation  (3.15): 


(4.2) 
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in  which  the  Pfl  dimensional  square  matrix  M'  has  elements; 


M’  - 
rs  ^ 


(4.3) 


We  taiow  from  prior  analysis  |A, 


(l/2N^)^ 


^2  I , (see  equations 


(3.11)  - (3.13)).  We  now  wish  to  obtain  rigorous  bounds  on  ®^id 


Aj ,A^ , . • . ,Ap^^ , in  the  limit  as  approaches  k 
is  used  to  Invert  (4.2): 


z Bin 


Cramer's  rule 


A^  - D^/Det!!M'l!  (4.4) 

in  which  D_  is  the  determinant  of  a matrix,  of  which  the  column  vector 
of  the  r.h.s.  of  (4.2)  is  the  rth  column,  while  all  other  columns  are 
the  same  as  in  M'.  Explicitly, 


u,(l) 

0 

0 

1 > • 

^2 

1 

1 * f 

r 

^1 

C2 

1 

1 

1 ‘'<•1 

,-i 

,-i 

r-1  1 

C2 

C2 

£3  - 

'1 

^3  1 

'> » 

• 

• 

# 

» 1 
• 

The  numerator  factors  into  Uj^(l)  times  a Vandermonde  determinant  of  the 

type  (3.17),  denoted  F, . In  the  limit  k — k , when  , the 

i z z min  i z 

denominator  (also  a Vandermonde  determinant)  vanishes  as  a linear  power: 


lim  DetllM'll  ■ (sin  k a)C.  (4.6) 

k — k ^ 

z z min 

as  the  first  two  columns  become  identical  and  » exp-ik^a,  with 
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Thus,  substituting 


k , “ 0 or  tt/s  . C,  is  a constant  near  k , . 

z min  1 7.  min 

the  knovm  value  of  yields; 


11m. 

k — k . 

z z min 


I sin  k^al  !c^/F^l (1/2N^) 


(4.7) 


The  vanishing  of  in  the  limit  is  ensured  by  the  non-vanishing  of 

the  Vandermonde  determinant  . 

In  the  limit  k — k . , the  analogous  expression  for  A„  differs 
z z min  B r 2 

from  (4.5)  only  by  interchanging  of  the  second  and  first  columns  in 
the  numer.’tor.  This  establishes  a second  result; 


lim.  A - -A  - K1/2N  y, 

k -k  ^ ^ ^ 

z z min 


(4.8) 


the  factor  of  iis  to  make  U real.  Finally,  for  r > 3 we  have  D^.  = 
Uj^(l)  X F^,  with  F^  appropriate  Vandermonde  determinant,  and  so: 


llm.lA^I  •=  (1/2N^)^1  F^/F^j  • (1/2N^)^  sin  k^e  - 0 (4.9) 


in  which  B is  constant  in  the  neighborhood  of  k , . We  have  made 
r z min 

use  of  the  fact  that  for  r > 3,  F = (C.”Co)  “ -21  sin  k a,  with  k — 0 

— r *^1  2 z z 

or  tt/b  . 


This  last  result  (4.9)  is  important  in  that  it  sharply  limits  the 
extent  to  which  the  evanescent  (decaying)  waves  contribute  to  the  solu- 
tion at  the  extremum.  It  establishes  the  qualitative  similarity  between 
all  the  Q“1  type  extremum  solutions,  for  arbitrary  band  structure,  with 
one  another  and  with  the  simplistic  P«1  limit,  equation  (4.1). 


(b)  Turning  to  the  more  complicated  analysis  of  Q > 2,  we  write, 
as  before,  ? "ZBC  equations"  for  m “ -(P-1),...,0  together  with  addi- 
tional equation  for  Uj^(ti^l).  We  wish  to  establish  that  the  functions 
U(m)  associated  with  the  extrema  (such  as  3nd  iri  Figure  3) 

vanish  at  the  extremum  (denoted  *)  while  other  functions  distant  from 

(3) 

the  extremum  but  degenerate  with  it  (such  as  point  in  Figure  3) 

have  finite  amplitude  (ni)  at  finite  m > 1. 

At  our  disposal  are  P-Q  distinct  complex  k's  and  Q+1  reals,  the 
latter  to  be  chosen  from  among  the  2Q  reals,  as  follows;  pick  one  pair 
straddling  a minimum  (such  as  ®nd  in  the  Figure),  and  choose 

the  remaining  Q-1  reals  all  unpaired,  from  among  the  positive  and  nega- 
tive k's  at  our  disposal.  The  unpaired  k's  may  be  lumped  together  with 
the  complex  and  imaginaries  for  the  purpose  of  an  analysis  such  as  in 
part  (a),  and  one  recovers  all  the  results  (4.2)-(4.9),  with  a slight 
-modification:  sin  k (h  - k ^^^)a  replaces  sin  k a on  the  r.h.s. 
of  Equations  (4.6),  (4.7)  and  (4.9).  At  an  absolute  maximum  or  minimum 
it  is  impossible  to  choose  Q-)-l  reals  without  including  at  least  one 
such  pair,  therefore  all  possible  Uj^(m)  vanish  at  finite  m as  the  energy 
approaches  its  extremum.  For  local  minima,  however,  such  as  the  case 
illustrated  here, it  is  possible  to  pick  a set  of  k's  such  that  no  pair 
straddles  the  minimur.,  all  remaining  distinct  as  the  energy  approaches 

the  relative  extremum.  For  example,  pick±k^^^^  and  of  which 

★ 

remain  distinct  as  e is  decreased  through  the  minimum  value  e . In  this 
case,  the  Vandermonde  determinants  all  remain  finite,  and  consequently 
the  amplitudes  and  U(m)  are  all  finite. 


Thus  solutions  straddling  a minimuni  (or  a maxiirun)  vanish  near  the 


surface  as  the  ninimuni  or  naxtmuni  is  approached,  whereas  other  solutions 
may  remain  finite. 

In  the  calculation  of  the  surface  l.d.o.s^®  as  in  Equation  (3.19) 
the  extrema  — local  as  well  as  absolute  — will  contribute  little  be- 


cause of  the  vanishing  of  the  amplitudes.  Two  interesting  consequences 

follow  for  the  quantities  p°  in  (3. 19):  first,  p°  at  the  band  edges 

...  . , min/maxi3/2  , , tnin/max,^ 

will  vanish  as  Ic-s  i rather  than  as  E-e  ; , the  latter 

n ' ' n ■ 

being  the  usual  behavior  near  the  band  edges  foi  the  bulk  edges  for  the 
bulk  d.o.s.  such  as  (3.20).  Second,  the  familiar  spikes  and  cusps  in 
the  d.o.s.  known  as  van  Hove  singularities,  caused  by  saddle  points  in 
the  energy  surface,  will  not  appear  in  the  1. d.o.s.  (3.19)  due  to  the 


vanishing  of  the  amplitudes  near  the  surface  at  these  stationary  points . 

We  have  proved  that  the  unperturbed  surface  1. d.o.s.  p"^  is  nar- 

n 

rower,  smoother,  and  (because  the  area  must  be  constant)  has  a higher 
maximum  at  the  center  of  the  band  than  its  bulk  counterpart.  As  m in- 
creases, the  structure  characteristic  of  the  bulk  progressively  re- 
appears . 


V.  EJ1AC7LY  SOLUBLE  SEMPLE  MODEL  PERrURBATION 
A perturbing  Hamiltonian  at  the  surface  will  mix  the  eigenstates 
U^j^(m)  of  the  unperturbed  ZBC  sem.i-inf inite  solid.  If  the  perturbation 
extends  to  P different  surface  planes  and  mixes  L distinct  energy 
bands,  the  solution  of  the  scattering  problem  will  require  the  diagon- 
alization  of  a (P  x L)-dimensional  matrix,  and  the  bound-state  ("sur- 
face-state") problem  requires  the  solution  of  an  equally  high-dimen- 
sional transcendental  secular  equation.  In  this  Section  we  solve  the 
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stnplest  non-trivial  ttodel  of  a perturbing  Hamiltonian  which  is  dia- 
gonal in  the  Wannier  representation  (L“  1)  and  confined  to  the 
first  -c“l  surface  plane  (hence  "P*  1).  It  goes  without  saying  that 
has  the  translational  sycanetry  of  the  lattice  in  the  x-y  plane;  the  com- 
plications which  result  from  breaking  this  synraetry,  as  in  the  example 
of  a single  ad-atom,  are  deferred  to  a subsequent  paper. 

The  matrix  elements  of  H in  the  terminated-solid  representation 

s 

are : 


^^s^nk.n'k' 


(5.1) 


in  which 


*Tk 


x,y- 


, i(k  X.  + k Y ) 

I e * J ^ J U , .(m,)0_(r-R.)  (5.2) 


nk  j n j 


'j 


and  Z.  ■ m a.  For  H diagonal  in  the  representation  of  the  Wannier 
J J ® 

functions  0 ("r-R.)  , and  localized  at  m=l,  the  matrix  element  is  indepen- 
n J 

dent  of  k and  k ' : 
z z 


,k  ''k  ,k  '‘n.r.’ 

’ x X y y 


(5.3) 


in  which 


8n 


- 


(r-Rj)H,(r)0„(r-Rj)6^ 


,1 


(5.4) 


The  eigenfunctions  of  the  total  Hamiltonian  -<•  , subject  to  ZBC, 

consist  of  scattering  states  and  possibly  bound  states,  also  known  as 
"surface  states".  We  start  by  constructing  the  latter.  Using  the  short- 
hand T “ (n,kn)  we  denote  these: 


25 


r 

I 


Schrodinger 's  equation  for  the  bound  state  is; 


(Hq  -t-  « Ej,  (5.6) 

We  make  use  of  the  knowledge  that  H , “ c (k  )V  , with  c^(k  ) a 

° o nk  t'  2 nk  t z 

•more  suggestive  -notation  for  the  band  energy  € (k)  . (We  also  occas- 

n 

M 

sionally  write  U.j.  (®)  to*  ^nk^°^ ''  Multiplying  Schrodinger ' s equation 
on  the  left  by  ®^d  integrating  over  electron  coordinates  yields 

a set  of  equations  for  '• 

<5-’> 

where 

A-  I U .,(l)F,(k'  ) 

. , nic  7 2 

z 

Self-consistency  requires  either  A • 0,  in  which  case  is  not  normal- 
izable and  thus  does  not  exist,  or  A V 0,  in  which  case  A be  fac- 

tored out  to  obtain: 


I - 


(5.8) 


where 


S^(E)  s I 


k-  e,(k^’)-E 


a 

- I-  f dk 
2n  J , 2 

— TT/a 


(5.9) 


Because  | u|  ^ , S^(E)  is  a quantity  0(1)  in  the  thermodynamic 

limit.  It  is  also,  in  general,  complex,  when  E overlaps  the  continuum 
of  between  the  relative  band  extrema.  Thus  for  (5.8)  to  be 


'1 
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satisfied,  must  lie  outside  each  continuum.  This  requirement  does 

•not  preclude  the  bound  state  at  one  value  of  t overleoping  the  con- 

tinuum  c . Ck  for  a different  t' . as  H does  not  mix  different  t's. 
T z ' s 


For  positive  (repulsive  perturbation)  E,  must  lie  above 

C_(k  ) E c'”**(k|i)  . But  (-S  (E) ) achieves  its  maximum  at  E •= 

T z max  n t 

)»  therefore  by  equations  (5.6)  there  can  be  no  surface  state 

T z max 

of  a repulsive  potential  unless  exceeds  the  critical  value 


(5.10) 


Estimates  of  the  numerator  of  S,,  from  either  (4.1)  or  (4.16)  and  ex- 
pansion of  the  denominator  of  S about  k using  e (k  + q)  a- 

e,(k  ) - Zlct^^^ia  ^ sin^(qa/2)  lead  to  the  above  estimate  of  the 
T'  2 max  zz  ' '.1 

T.h.s.  of  (5.10).  (See  Appendix  A.)  Evidently  the  instability  against 
bound-states  is  greatest  at  points  in  the  2DBZ  where  the  curvature  at 
the  relative-band  maximum  is  the  smallest  (i.e.  at  points  t where  the 
effective-mass  is  the  greatest). 

For  negative  (attractive  perturbation)  the  bound  surface  state 

energy  E lies  below  e (k  . ) k e™^(k||)  . We  find  there  can  be  no 
T T 2 min  n H 

such  state  unless  the  magnitude  of  the  attractive  potential  exceeds  a 
critical  value; 


(5.11) 


with  the  estimate  of  the  r.h.s.  being  obtained  in  a similar  way  to 
(5.10).  (Note  that  S is  negative  in  the  above.)  Where  • 0 any 

arbitrary  weak  perturbation  will  cause  a surface  state  to  '‘peel  off"; 


see  Figures  1,  2 and  4, 
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If  a bound  state  exists,  it  must  be  normalized.  One  readily  finds 


for  the  normalized  amplitudes; 


k(1>  i 


(5.12) 


in  which  the  derivative  is  to  be  taken  using  (5.8). 

Note  that  may  be  sufficiently  large  to  cause  surface  states  for 

some  k|j  but  not  for  others.  While  in  the  n.n.s.c.  (lOO)-band-struc- 

ture  equations  (5.9)  and  (5.10)  are  inaens itlve  to  kjj , in  the  (Oil) 

orientation  of  the  same  band  structure,  the  same  a at  k a r:  ± r:2  ^ 

22  y o 

results  in  instabilities  shown  in  Figure  4. 

Scattering  theory  gives  us  the  form  of  the  remaining  states; 


K u - Y ^(r)  + I 6^  , .8  , . L,  ..Y  , ,(r) 


T,k  'nk 


utj.  k ,k  ' k ,k  ‘'k,k'‘nk' 

kVk  3c’  X y y 


(5.13) 


As  in  the  analysis  of  the  bound  state,  we  write  Schrodinger 's  equation 
for  this  eigenftinction : 


<«0*V*T,fc 

’2  2 


(5.14) 


or,  more  explicitly: 


e(k  )Y  -(r)+  I y 

*^2  k 'Vk  z’  2 2 

2 2 


2 *^2  z 2 2 2 2 2 2 


E(k)r,  <r),  r L,  .(Dl 

2 *^2  ^ ^ * 
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and  equate  coefficients  of  the  Yj^’s  (the  subscripted  T and  kj|  omitted  for 

typographical  clarity).  Starting  with  k we  find: 

z 


E^(k^) 


'Tk 


(1)1 


'Tk 


k,'Vk.  z z z 


(5.16) 


writing  for  For  most  purposes  E (k^)  can  be  replaced  by  e(k^) 

which  it  interlaces.  However,  if  it  is  desired  to  calculate  the  cor- 
rection, the  above  may  be  used  to  do  so  with  the  aid  of  the  cofficients 
k'"  '^®  obtain  by  equating  coefficients  of  'j'k^ ' (*^2  ' ^ H)  • 


^^V^k.k' 


Sn^k< 

z z z 


vhere 


^k 


z 


These  are  readily  solved,  to  obtain: 


(5.17) 


'k  ,k  ' 
z z 


* y-'**  8„5,(E,(y>T' 


E,(k,)  - .,0^, 


(5.18) 


For  most  applications  E.,(k^)  can  be  replaced  by  Its  thermodynamic 
limiting  value  £.^(k^)  in  this  formula.  Note  that  the  denominator 
(1+  g„S)  cannot  vanish,  as  S is  complex  in  the  rsnge  of  interest  owing  to 
the  branch  cut  along  the  continuum  of  the  relative  band.  However,  if 
there  is  a bound  state  outside  the  continuum  then  at  a symmetric  point 
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with 


= In 


%(Rf2)  - [i(R+2)^-n^ 


2 , tnax/min.  -1 


and 


^“zz  ^ '^T~®T^*^z  max/mir?' 


(5.23) 


R being  effectively  the  ratio  of  binding  energy  to  band  width  for  the 
bound  state.  This  estimate  does  not  replace  a quantitative  calculation 
based  on  actual  band  structures,  but  can  serve  as  a guide. 

We  can  now  compute  the  surface -state  contributions  to  the  local 
density-of-states  (l.d.o.s.),  at  fixed  t » (n,k|i): 


PT.m;bd.st.f")  “ 


(5.24) 


Proceeding  with  the  scattering  states  in  the  same  way,  we  write: 

i(k  X.  k Y.) 


with 


-i 

i,  , - N ^ I e 

Jt.y  R. 

^ j 


* ^ ^ ^ W,  u (n,)0  (r-R.) 

j'  n"  3 


(5.25) 


“T.k  <">  • “T.k  <">  * H ,k  '”T.k  ■<")  (=■“> 

z z z z ’ z 


according  to  (5.13).  Using  given  in  (5.16)  we  obtain: 


*n”T.kW  "T,k  ''''"T.k  , W 

Wt.k^<">  ■ ”.,k/”)  - i;i;;5;77r  kt'  ~ 


• * • 


(5.27) 
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abbreviating  c^(k  ) by  e.  This  simplifies  at  nt“l: 


“t.k  <»  ■ "T.k 


(5.26) 


Again  the  m-dependent  integral  can  be  estimated  for  ic  > 1 (see  Appen- 
dix A)  resulting  in  the  following  estimate  of  (5.27): 


W ^ (m)  • (2/N  sin  ((k  -k  . , )ma  - 6,  ) (5.29) 

z'  '^2  zmin/max^  k^ 


The  scattering-state  (i.e.  bulkr)  contributions  to  the  l.d.o.s. 

are; 


^.m.sc.st.^-^  ' I l''T,k  6(E-c,(k^))  (5.30) 


(5.24)  and  (5.30)  have  been  obtained  at  fixed  t - (n,k).  The  results 
are  sketched  in  Figure  4b.  The  following  sum  rule  follows  from  complete- 
ness ; 


1 - r dElp  . (Z)  -f  p (E) 

J '^T,m,bd.st  T,n,8c.st.  ' 


(5.31) 


We  also  define  the  total  band  l.d.o.s.  by  the  sum  over  k„ 


^n,m,bd.st.^^^  ^x,y  ^T,m,bd  .s  t . 


and 


-1 


^n ,tn , • c ..s ^*.y  k,|  .st . 


(E) 


(5.32) 
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as  sketched  in  Figure  4.  These  p's  also  satisfy  the  sum  rule  (5.31). 


In  the  limit  g “0  the  bound-state  d.o.s.  all  vanish  and  the  scattering- 
n 

state  d.o.s.  become  equal  to  the  unperturbed  functions  defined  in  (3.19) 
snd  (3.20),  both  of  which  have  unit  area.  In  the  opposite  limit  (large 
Igjjl)  it  is  possible  for  the  scattering  contributions  to  become  small 
and  the  area  sum  rule  to  be  saturated  by  the  bound-state  contributions 
^ this  is  the  Interesting  lim.it  of  complete  "surface  bands". 

VI.  VACUUM  STATES 

For  the  sake  of  completeness,  the  set  of  eigenstates  for  the  solid 
must  be  augmented  by  a set  of  states  for  the  outside,  the  'Vacuum 
states".  So  far,  we  have  used  only  the  Wannier  functions  of  atoms 
within  the  solid  as  our  basis.  Let  us  now  test  their  completeness; 

Z fl*(r'-R.)0^(r-R.)  • A,  (r'.r)  (6.1) 

Rj(m>0)  « J J + 

all  n 

Deep  in  the  bulk  (z  or  z'  > Pa,  with  P the  cut-off  integer  introduced 
in  a preceeding  section)  we  have  the  limiting  behavior: 

lim  L (r',r)  — £(r'-r)  (6.2) 

z.z'— * 

i.e.,  the  set  of  Wannier  functions  of  the  solid  become  complete  for 
the  description  of  arbitrary  functions  deep  in  the  bulk.  But  near  the 
surface,  the  function  A_^(r',r)  differs  — perhaps  seriously  so,  for 
some  applications  --  from  the  required  delta  function.  Obviously, 
vacuum  states  are  required  for  the  description  of  electrons  outside  the 
solid.  Less  obviously,  they  are  also  required  to  augment  the  states 
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of  the  solid  for  a good  description  of  the  surface  region,  both  inter- 
-nal  and  external  to  the  solid. 

Evidently,  one  can  re-introduce  the  missing  functions  In  a manner 
analogous  to  those  which  were  retained,  and  define: 

r a*(r'-R  )0  (r-R  ) e L (r'.r)  (6.3) 

R . (m  ,<0)  " J " J 

J J- 

ali  n 


Now , 


i^(r',r)  * i_(r',r)  » i(r'-r)  (6.4) 

everywhere . But  as  the  contribution  of  t becomes  vanishingly  small 
when  r ' or  r are  deep  in  the  bulk,  so  does  the  contribution  of  be- 
come negligible  at  a comparable  distance  outside  the  solid.  It  is  thus 
mainly  within  the  surface  region  that  both  sets  of  states  must  be  used, 
and  we  see  that  the  surface  region  possibly  has  greater  structure  then 
might  have  been  anticipated;  see  Figure  5. 

Well  outside  the  solid,  the  are  badly  suited  for  describing 

free  electron  motion.  The  perturbances  due  to  the  presence  of  the 
solid  do  not  manifest  themselves  until  the  surface  region  is  reached. 
Thus  for  the  vacuum  states,  plane  waves  are  required,  at  least  asymp- 
totically. We  wish  to  construct  plane  waves  using  the  0 's,  so  as  to 

n 

ensure  the  orthogonality  to  the  states  inside  the  solid.  The  simplest 
procedure  consists  of  defining  the  following  basis  functions: 

Cj^(r)  c fd^r'  £_(r',r)  (6.5) 
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in  which  the  normalization  constant  is: 


C, 

K 


^ik-(r--r-)^  (r-.r") 


(6.6a) 


In  the  limit  that  the  volume  of  the  vacuum  region  D k (L  N a N a ) is 

2 X X y y 

Infinitely  greater  than  the  volume  of  the  surface  region  0(Pa  N a N a ) 

2 X X y y 

most  of  the  contribution  to  this  integral  comes  from  the  asjicptotic  re- 
gion where  t 6{r'-r")  and  therefore. 


C ,«  L N a N a - T., 
k 2 X X y y 


(6 . 6b) 


and  is  independent  of  k. 

Variationally , the  energy  of  each  state  is: 


(6.7a) 


Because  any  perturbations  Hg  are  limited  to  the  neighborhood  of  the  sur- 
face, and  because  is  essentially  a plane  wave  (except  in  that  same 
■neighborhood)  we  have; 


Ej^  - ^^k^/2n+  0(Pa^/Lp  (6.7b) 

Although  normalized,  orthogonal  to  the  states  of  the  solid,  and 
practically  "sharp"  in  energj’,  the  states  suffer  from  two  Important 
defects.  The  first  concerns  their  lack  of  mutual  orthogonality.  This 
is  serious,  for  it  would  prevent  us  from  remedying  the  second:  the 

fact  that  they  are  not  exact  eigenstates  of  the  total  Hamiltonian,  i.e., 
that  a perturbation  Hg,  existing  near  the  surface,  may  mix  them  with  one 
another  and  with  the  states  of  the,  solid.  The  application  of  scattering 
theory,  the  calculation  of  matrix  elements  such  as  are  required  in  the 
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photoelectric  effect  and  other  surface-related  dynamical  phenomena , 
problems  in  chemisorption,  etc.,  all  require  the  use  of  an  orthonormal 


basis . 

For  this  purpose , Schmidt's  orthoeonalizatlon  procedure  fails. 

This  becomes  izmediately  apparent  as  one  attempts  to  use  it.  One  starts 
vith  a function,  say  orthogonalizes  the  remaining  functions  to  it 

normalizes  the  new  singles  out  the  next 

function  — say  — and  repeats  the  process,  iterating  until  the  set 
of  functions  is  exhausted.  As  this  procedure  is  in  the  nature  of  an 
algorithm,  it  would  appear  ideal  for  computer  applications;  however, 
suppose  we  wish  to  determine  the  new,  properly  orthonormalized , func- 
tion replacing  C at  a finite  energy  Z above  the  bottom  of  the  con- 

K K 

tinuum.  This  function  can  only  be  reached  after  an  infinite  number  of 
steps!  (This  is  not  an  unimportant  objection,  because  it  is  not 
remedied  by  box-quantization,  although  this  reduces  the  number  of  steps 
to  a large,  but  finite  number  of  steps.  It  is  not  at  all  clear  that 
the  orthonormal  functions  approach  limit  functions  as  the  size  of  the 
box  is  allowed  to  approach  infinity.)  One  desires  a formulation  capable 
of  yielding  the  states  at  an  energy  Z^  without  the  need  for  calculating 
all  the  states  belonging  to  Z,  ' < E,  first.  A second  objection  is  that 

K 

the  Schmidt  procedure  is  totally  arbitrary.  If,  instead  of  starting 
with  we  had  started  at  some  arbitrary  the  function  which 
this  procedure  ultimately  would  yield,  would  differ  greatly  from  the 
initial  and  indeed,  eve ry  function  obtained  in  this  way  would  dif- 
fer from  those  obtained  by  the  initial  prescription. 
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For  our  applications,  we  have  therefore  developed  a procedure 
■which  is  specifically  tailored  to  contiauun  states.^  The  Schmidt  pro- 
cedure may  continue  to  be  used  for  any  finite  subset  of,  e.g.,  local- 
ized states,  while  the  procedure  we  now  outline  will  be  both  practical 
and  well  suited  for  the  continuum  states.  The  arbitrariness  in  the  re- 
sults we  have  noted  will  be  parametrized  by  a set  of  quantities  c^, 
which  we  can  choose  arbitrarily. 

We  first  -make  use  of  the  obvious  identities: 

C(r\r)  - i_(r,r')  and  Td  r"  A*  (r  ,r")i^  (r  ' , r")  - ^,(r',T)  (6.8) 


together  with  (6.4)  to  obtain  the  overlaps  in  the  form: 


Jv  “ \.k' 


- (1/L  ) ZT  (k')F  (k) 

all  n 


with: 


F^  _ (k)  E a ,/n)^  fd,r  0„(r-R.), 

IlfRt  Z I)  j 

j 


(6.10) 


obtaining  (6.5)  in  the  form: 


C.  (r)  - (1/H)^  - (n/l,)^  Z T*  _ (k)0  (r-R  )!  (6.11) 

^ Rj(m>0)  ri.Rj  n j 


The  vacuum  region  outside  the  solid  is  the  given  volume  ■ 

L N a N a , with  L ultimately  allowed  to  become  Infinite.  Thus  any 
zxxyy’z 

finite  length  such  as  Pa  may  be  neglected  bv  comparison.  The  arguments 

z 
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given  following  equations  (6.5)  and  (6.6)  indicated  that  the  wave  func- 


tion (6.11)  vanishes  for  z > Pa^,  thus  only  tenr.s  with  < P need  be 


retained  for  L < z < Pa . 
z 


We  now  start  our  procedure,  by  a simple  step,  in  which  only  a 
single  function  <2  (r)  is  involved. 


Instead  of  (6.11)  let  us  first  assume  that  the  basis  set  consists 


.(0) 


of  C,  defined  as  follows: 


tWhr)  . a/n’>  1.“'"  - f'“>* 


(6.12) 


in  which  the  0(r)  to  which  all  are  orthogonal  is  nonr.aliced  and 


E ^e  C (r)d.,T,  which  is  0(1). 

k >>  0 J 


The  following  results  tnay  be  proved  using  completeness  of  the  e 


(6.13) 

ik.r 


(0),2 


(1/n)  I IF^  = 1 


(i/r.)  I - 0 (r) 

k k O 


(1/n)  z F^°^  “ 0 


(6.14) 


Our  procedure  is  to  construct  a pseudo-Kamiltonian  in  the  space  of  the 
functions  (6.14)  having  0(r)  as  its  ground  state;  the  excited  states 
will  then  be  the  desired  orthonormal  set.  So  let  us  define  the  matrix 


■elements  of  H • plane  wave  representation: 


V . Ol  fi  - i r*r 
1^1.  1.1  <1—  ° 1,  1, ' ' 


(6.15) 


‘k,k’  2m  k,k'  n k'^k' 
with  \ an  adjustable  parameter  and  to  be  determined.  Evidently,  the 


1 
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bound  state  (ground  state)  takes  the  form 


* - (l/n)  Z g,  e’ 


(6.15) 


■nd  has  energy  - c . Schrodinger 's  equation  for  the  and  leads 


/tiid  h.  r*  T-  r 

^ 2m  ^o^®k  “ n ^k  ^i^k’^k' 


(6.17) 


The  non-trivial  solution  of  these  equations  requires; 


(1/n)  Z 


“t ^ 

2m  o 


(6.18) 


and  thus  the  normalized  solution  is : 


’k  2 2 


V ,2 

2m  ^o  2m  ^o 


(6.19) 


Comparison  of  (6.16)  and  (6.14)  indicates  that  we  need  8j^  “ , the 

overlap  function  defined  in  (6.13).  Thus,  we  choose: 


c;;-  2 1.2/20+ 


(6.20) 


which,  together  with  (6.18)  determines  H'.  Note  that  remains  arbi- 
trary ; it  in  fact  represents  the  extent  of  arbitrariness  in  the  pro- 
cedure, and  parametrizes  the  infinite  number  of  possible  orthonormal 
sets  which  can  be  constructed  having  the  desired  properties.  Using  any 

convenient  criterion  to  determine  t > 0,  or  arbitrarily  setting  it 

o ” 
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zero,  we  can  construct  the  corresponding  "perturbation  pseudo-Katr.ilton - 
ian"; 


r 

I 


K 1,1 

k,k 


1 + e^)(Fj[°^(A^k’^/2i:i+  e^) 

^ (1/r.)  I + Eq) 


(6.21) 


The  next  step  is  to  constrict  the  scattering  states  of  H*.  They 

ere : 


.,(r) 


where 


\ir) 


i ik*r 

ie 


■+ 


(1/n)  I 

k' 


ik’-r 


1 


(6.22) 


We  have  already  solved  the  scattering  probleiE  for  a separable  H'  in  the 
preceeding  section.  The  xesults  are; 


.(0)*,2,  2 


•(0) 


(ft"k^/2!n+  e (^l^k^/2ir, 

O K 


2m 


(k^  - k'2) 


^(0) 


(6.23) 


where  Dj^  is  complex: 


e ) ^-  - 

'■  2m  o'*  n k"  ' ^k"  i ^2 


^^k"^/2m  + e 


(6.24) 


This  specifies  the  orthonormal  set  of  wave  functions  (6.22),  positive- 
energy  eigenfunctions  of  the  synthetic  Hamiltonian  (6.15). 

To  further  orthogonalize  this  infinite  set  of  functions  to  a given 
flj^(r),  we  merely  iterate  the  procedure,  replacing  exp(ik-r)  in 
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(6.12)  - (6.16)  by  *j^(r),  i.e. 

Cj[’\r)  E (1/n)^  i»j^(r)  - 0^(r)!  (6.25) 

■nd 

- (l/.Q)E^L^°j^,  e"^'^’’’'!0^(r)d3r  (6.26) 

differs  from  (6.10),  it  should  be  noted.  If,  «s  we  shall  assume  0^^  is 

13 

orthogonal  to  0 and  is  normalized  as  well,  ^ then: 
o 

a/njElF™!"  • 5 i JS'S'' 

•=  Jd3rd3r'[^^  I *^(r'):^^(r)  + 0*  (r ')  0^  (r)  10*(r)  0^  (r ' ) 

- Jd^rd^r'  6 (r-r ')0* (r) 0^ (r ') 

- 1 (6.27) 

in  which  we  have  added  a term  which  vanishes  (due  to  the  orthogonality 
of  0^  and  0^^)  and  used  completeness.  This  proves  the  equivalent  of  the 

first  line  of  (6.14).  To  prove  the  analogue  of  the  second  line,  we  also 

add  a zero  term  containing  0^  and  use  completeness,  obtaining: 

(1/P.)  Z f5^N.  (r)  - 0 (r)  (6.28) 

k K K i 

«nd  finally,  orthogonality  of  the  new  functions  to  the  0^  ensures  that: 

(l/u)  I ■=  0 (6.29) 
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The  TDfltrtx  elements  of  the  new  pseudo-Hamlltonlsn  In  the  representation 
of  the  are; 


\,k'  “ '^k*k,k'  ' n SS' 


(6.30) 


where  * #l^k^/2m  + 0(l/f!)  , and  X and  Gj^  for  the  separable  Interaction 
must  be  deterrcined , as  previously,  so  as  to  make  0^^  the  unique  bound 
state,  and  the  remaining  states  Ipso-facto  the  orthonormal  set  of  posi- 
tive energy  states  orthogonal  to  0^^.  By  strict  analogy,  one  obtains 
that  these  states  take  the  form: 


for  -L  < z < Pa  , 
k z z 


where 


♦ ">(r)  . (1/a)  I 1<”,  *^,(r) 


(6.31) 


In  which 


(1)  _ ^k  ^ ^^k  ^k'^  ^^k’  ]_ 

H,k-  - c,  - t,.  • ,(1) 

k 


(6.32) 


where 


(c,  + e.)  r- I.  F.„  7 7 

k 1 n k"  ^ Cck"  " ‘^k^ 


(6.33) 


In  the  Iterative  scheme,  exp(ik-r)  played  the  role  of 
(r) , etc.  We  now  Iterate  any  finite  number  of  steps,  and  re-express 
the  results  in  terms  of  the  original  plane  waves: 


(1/0  I,  e^‘^’- 


(6.34) 


.Ji 


In  which  the  functions,  form  the  desired  orthononnal  set 

orthogonal  to  the  set  of  0 ,0,  The  recursion  relations  are: 

o 1 n “1 


14 


and 


w(n)  (n-1)  „(n-l)  i _ . (n-1)  (n-1) 

\vc’  Hk’  ^ \k«  ^ n 1..  \k"  \"k* 


(6.35) 


.(n) 

^ 


'k  ■ 'k' 


,(n) 


(6.36) 


®n^  C k"  ' (£k"  “ 


>(n) 


(6.37) 


(6.38) 


The  € 's  are  arbitrary  but  should  be  chosen  O(work  function)  for  real- 
n 

2 7 

Isa.  The  "energies"  are  Cj^  * k /2m.  While  for  the  nonr.ali- 

ration  n , the  energies  and  the  procedure  outlined  above  all  have 
to  be  modified  for  practical  purposes,  we  orthogonalize  to  a finite 
number  of  planes  P in  the  solid,  belonging  to  L bands,  so  n is  in  fact 
finite  and  < L x P. 


VII.  MORE  GENERAL  SURPACE  PERTURBATIONS 
The  model  perturbation  treated  in  Section  V was  exactly  soluble 
but  hardly  realistic;  effects  of  a surface  can  be  expected  to  extend 
well  beyond  the  first  surface  plane.  Although  in  scattering  theory, 
few  problems  are  soluble  in  closed  form,^^  we  have  fortunately  found 
two  forms  of  surface  perturbations,  both  extending  to  P surface  layers 
and  involving  a number  L of  different  bands,  for  which  a complete  and 
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and  explicit  solution  can  be  constructed  in  a finite  number  P x L steps. 
They  are : 

(a)  diagonal  in  Wannier  function  representation,  l.e.  having 

band "diagonal  matrix  elements: 


^«s^k,n’k'  “ 8n^“>«n.n-\,,.k'„Unk<“'^'^n-k-(“> 


(7.1) 


in  the  terainated-solid  (ZBC)  representation.  This  H may  be  inter- 

s 

preted  as  merely  a local  perturbing  potential  acting  separately  on 
«ach  plane  and  for  each  band  energy.  It  nevertheless  represents 
a formidable  problem,  the  solution  of  which  requires  special  methods 
outlined  below. 

If  the  perturbing  potential  varies  rapidly  across  an  individual 
surface  plane,  it  perforce  has  interband  matrix  elements.  This  feature 
is  lacking  in  model  (a) . But  introducing  it  can  easily  involve  us  in 
huge  matrix  diagona lira Cion  which  we  wish  to  avoid.  Still,  one  would 
like  to  be  able  to  study  the  effects  of  interband  transitions  on  such 
important  properties  as  recombination  of  electrons  and  holes  at  the  sur- 
face on  semiconductors,  optical  selection  rules  at  the  surface,  etc. 
Therefore  we  develop  a second  model  in  which: 

(b)  mixes  different  bands  but  in  a "separable"  manner.  The 
matrix  elements  in  the  ZBC  representation  are: 


(7.2) 


The  band  index  n may  also  include  the  positive-energy  (vacuum)  states, 
provided  a suitable  high-energy  cut-off  is  introduced  compatible  with 
the  maximum  slope  of  across  a plane. 
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Although  the  two  soluble  models  given  above  represent  entirely 

different  physical  situations,  they  may  be  solved  by  a conanon  proced- 

2 

ure.  First,  note  that  \ (tn)  corresponds  to  g (m)  . Thus  for  repulsive 

n n 

potentials  we  take  v (m)  real,  for  attractive  ones  pure  imaginary.  It 
n 

* 

Is  also  convenient  to  define  an  operation  i « f for  all  complex 

•it  ^ ^ 

functions  except  v (m) , for  which  v “ V if  V is  real,  but  v * -v  if 
n 

V Is  Imaginary,  those  being  the  only  possibilities. 

We  now  show  that  (a)  is  just  a special  case  of  (b) . For  with 
2 

g (m)  * Y (m)  and  n limited  to  a single  value,  equations  (7.2)  and 
n n ~ 

(7.1)  coincide.  By  solving  for  each  band  ri  individually,  (b)  may  be 
used  to  solve  (a).  Thus,  it  suffices  to  study  (b) , the  more  general  of 
the  two  models . 

Our  procedure  now  is  related  to  that  of  the  preceedlng  Section. 

Instead  of  diagonalizing  H + T H (m)  all  at  once,  we  first  diagonalize 

o ms 

+ H^(l),  re-express  the  remaining  Hg(m)  in  the  new  eigenfunctions, 
diagonalize  (H^ -(■  Hg(l))  + Hg(2),  repeat  the  performance,  etc.  After 
P steps  the  total  Hamiltonian  is  diagonal,  and  all  the  eigenfunctions 
and  eigenvalues  — both  continuous  and  discrete  — are  known  in  the 
presence  of  the  complete  perturbation. 

Let  us  work  out  the  mth  step.  We  suppose  that  “ 


[H 


m-1 


^ ^ ^ H^(m')l  is  already  diagonalized,  and  now  wish  to  work  out 


the  effects  of  the  mth  perturbation,  H^ (m) . In  the  representation  of 
the  eigenstates  of  the  perturbation  has  the  matrix  structure: 


o'  “ t (m)t  , (m)6  . 

S q.q  q q 


(7.3) 
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vhere 


t (m)  * V X™)  (7.^) 

q q q 

and  q 8 quantuo  number  running  over  the  continuum  states  (n,k)  as  well 
and  the  bound  states  of  is  the  amplitude  of  the  cor- 

q 

responding  eigenstates  on  the  pth  plane.  Evidently, 


t (m) 

q 


0(N  S for  q in  the  continuum  of 
z 


0(1)  for  c ranging  over  the  bound  states 
of  ■ 


(7.5) 


The  eigenstates  of  the  new  Hamiltonian  h/'^^  e -f  (m)  are  ob- 

tained by  the  standard  methods  used  in  previous  Sections  V and  VI,  and 
result  in  new  amplitudes  given  by; 


q 


(p) 


M 


<p) 


q.q  q 


(7.6) 


with 


, (m) 
q.q 


t 


T , (m)  t (m) 


ri  + 


n q 


-1 


(7.7) 


and 


q q 


(7.8) 


The  normalization  constant  may  be  guessed  to  be  trivally  1 + 0(N^  ) for 
the  scattering  (continuum)  states;  in  general  It  is: 
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M - 1 + t (m)?"  (m)  ^ ln[;i  (E  )] 
q q q oEq  o q 


(w) 


(7.9) 


•■nd  the  new  energies  E 6 e ere  related  to  the  previous  eigenvalues 
(0-1)  . 

e.  « 'by: 

q q 


q q q q ® q 


(7.10) 


Thus,  by  (7.5)  for  continuum  states  the  energy  corrections  are  negli- 
gible, and  the  new  eigenvalues  interlace  the  old.  For  bound  states  of 
j^(o  1)  amplitudes  t^  are  0(1)  and  the  energies  are  shlfed  by  a fi- 
nite fraction,  and  equation  (7.10)  is  then  a non-trivial  transcendental 
equation  for  the  new  energies.  Finally,  new  bound  states  can  appear  as 
a consequence  of  the  perturbation  at  m.  These  will  appear  as  new  roots 
E of: 


1 (E)  » 0 


(7.11) 


at  energies  E which  differ  from  any  of  the  set  of  E 's  we  have  obtained 

q 

(either  bound-  or  continuum  states)  using  (7.10) . The  formulas  (7.6)  - 
(7.9)  for  the  eigenstates  can  still  be  used  for  the  new  bound  states; 
replacing  E^  by  E the  solution  of  (7.11)  in  these  formulas,  we  observe 
that  the  new  pole  in  L dominates  the  rest,  that  it  appears  both  in  the 
numerator  and  in  the  normalization  denominator,  that  the  terms 
t^(l+j/)  ^ will  cancel  thereby,  thus  leaving  a remainder  in  precisely 
the  form  of  the  normalized  bound  state  derived  in  (5.5)  and  (5.12). 

We  proceed  to  the  next  stage.  Setting, 
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(7.12) 


t (nn-1)  * V (nH-l)W^°^\n>f  1) 

^ q q 


•where  W^™\n>+-1)  the  amplitude  (7.6)  at  p • m^l,  one  solves  for  the 

q 

effects  of  the  perturbation  at  m -r  1.  Equation  (7.3),  with  m + 1 


■replacing  m,  -now  defines  Hg(m*-1).  The  procedure  is  repeated.  P such 

Iterations  following  the  initial  condition  W ^°\p)  = L’  , (p)  suffice 

q nK 

to  obtain  explicitly  the  exact  eigenstates  and  eigenvalue  spectrum  of 

(P) 

the  total  perturbed  Hamiltonian  K = H -f  H , and  the  corresponding 

os  r a 

CP) 

amplitudes  (p)  on  the  pth  plane,  for  any  1 < P S • 


VIII.  CONCLUSION 

In  this  paper  we  have  related  the  band  structure  to  the  notion  of 
the  chemical  bond,  have  shown  how  the  breaking  of  bonds  at  a surface 
alters  the  eigenstates,  and  how  additional  perturbations  can  create 
bound  states  among  these  — - identified  as  the  surface  bands.  The  ex- 
tinction of  the  bulk  van  Hove  singularities  in  the  surface  density  of 
states  is  an  interesting  new  result,  which  we  proved  rigorously  for  the 
terminated  solid.  The  introduction  of  ordinary  perturbing  Hamiltonians 
affects  surface  amplitudes  smoothly,  so  the  disappearance  of  van  Hove 
singularities  is  undoubtedly  a completely  general  feature.  Of  course, 
new  van  Hove  singularities  in  cases  where  there  exist  2D  band  of  sur- 
face states  can  be  expected,  but  they  are  then  surface-related  proper- 
ties. Our  theory  can  have  consequences  in  the  interpretation  of 
photoelectric  emission.  This  Important  probe  is  one  of  several  which 
are  specific  to  surfaces  and  require  a good  understanding  of  the  posi- 
tive-energy (vacuum)  states.  Fox'  the  latter  we  have  developed  an 
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orthonormalization  program  which  is  effective  for  continuum  states  and 
complements  Schmidt's  method  for  orthonormalization  of  discrete  states. 
We  also  recapitulated  our  earlier  criterion^  which  differs  from  those 
which  are  presently  well  known^^  for  the  instability  against  the  for- 
mation of  surface  states  against  arbitrary  perturbations.  Now,  let  us 
turn  to  the  nature  of  these  perturbations. 

Surface  perturbations  have  no  universal  character  and  so  must  be 
decided  on  a case-by-case  approach.  The  study  of  metals  will  differ 
from  Insulators,  and  both  are  different  from  the  large  variety  of  semi- 
conductor terminations.  Let  us  start  with  the  insulators;  we  have 
seen  that  the  eigenstates  are  in  one-to-one  correspondence  with  the  un- 
perturbed states,  so  that  if  we  populate  the  valence  band,  including 
surface  states  associated  with  it,  charge  neutrality  is  ensured.  How- 
ever, the  introduction  of  the  surface  has  invariably  contributed  to  a 
surface  dipole  moment.  (Electrons  have  been  removed,  the  wavefunctions 
of  which  extended  into  the  solid,  whereas  the  compensating  ionic  charges 
lay  outside.)  The  potential  associated  with  this  dipole  layer  repre- 
sents the  variation  in  the  Madelung  potential  — the  electrostatic  po- 
tential felt  by  the  electrons  at  any  site  — at  the  neighborhood  of  the 
surface.  The  potentials  and  the  charge  imbalances  they  engender  must 
be  related.  Poisson's  equation  and  the  Hartree  apprcxLmetlon  cay  a- 
chieve  this  relation,  as  may  the  Hartree-Fock  method.  Such  methods 
have  indeed  been  applied  to  surfaces,  and  a vast  literature  now  exists 
on  them.^^  Nevertheless,  applications  using  the  present  methodology 
may  have  advantages  over  the  previous  ones,  many  of  whlch^^  require  the 
ab-initio  solution  of  Schrodinger 's  differential  equation  and  make 
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little  or  no  use  of  the  information  already  available  from  the  study 
of  the  bulk. 

The  problems  of  metals  are  different  but  equally  challenging.  The 
■electron  states  are  popul ated  up  to  the  Fermi  energy , which  is  de termined 
by  the  bulk.  But  the  states  lying  below  may  have  small  amplitudes 
at  the  surface,  and  thus  the  surface  electronic  charge  may  be  deficient. 
We  see  that  this  would  result  in  an  attractive  corrective  potential, 
which  in  turn  would  increase  the  amplitude  of  the  states  at  the  sur- 
face, until  charge  neutrality  is  achieved.  Either  Hartree  or  Hartree- 

Fock  methods  may  be  used  here  or  the  more  modem  schemes  developed 

18 

by  Kohn  and  his  collaborators.  However,  one  suspects  that,  as  in 

19 

metals,  a judicious  application  of  Friedel's  sum  rule  may  prove  useful , 
decisive,  or  at  least,  a s implif ication . Ve  are  presently  investigat- 
ing this  possibility  as  well  as  those  several  other  issues:  catalysis, 

■ad-atoms,  etc.,  which  are  logically  related  to  the  present  work. 
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APPENDIX  A 


Near  an  excreisur;  U ~ sin  qa  and  e(k  ) p;  c , , * 2|o  I a"'^6ln^'ioa  . 

z max/nin  zz  ‘ ■ 

vhere  g e k - k . . We  use  these  expansions  not  just  for  scall 

z znax  rr.tn 

q,  but  throughout  the  EZ  to  estimate  such  quantities  as: 


n'^  I I*(l)U(ir)rE  - e(k^)l"^ 


(A.l) 


Thus,  the  above  is  approxinately : 

a - 11-22 

I • — ' dT|  sin  Tla  sin  Tlam:  A -i-  2 o 'a  sin  ?"a' 
ti  TT  >.  ' ' - zz 

— rr/a 


(A. 2) 


where,  for  E outside  the  continuum  (bound  state)  we  derive  (5.22)  as 
follows : let 


E - e. 


'max/min 


(A. 3) 


be  the  binding  energy,  and  R the  dimensionless  quantity 

R = 2 A a /]  O'  ].  With  J defined  by  I = Za'^J  /]  al  , we  have; 
zz ' m m m • zz 

I »rr  2-1 

~ j dT]  s in  T)  6 in  Tim  [ R + -4  s in  %7|  1 


- c”^^[(Rf2):-r-i]‘^ 


(A. 4) 


where  the  contour  is  the  unit  circle.  Evaluating  the  residue  at: 


%(Rf2)  - [i(R+2)^-ll^ 


(A.  5) 


«e  obtain; 


J - J,  exp-(m-l)v 

Q 1 


(A.  6) 
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vhere  is  expressible  precisely  in  terms  of  S^(E)  , and 
in  (5,23). 

For  E in  the  continuum,  we  have  replaced  E by  e (k^) 


J •=  J e 
tn  o 


iqam 


-where  J is  real  and  J,  « S (e(k  )).  Thus, 
o 1 ' 2 


S (e  (k  ) ) a:  1 S,  exp  ica 

T Z I 


For  the  derivation  of  (5,29),  we  use 


U ss  (2/N  )^  sin  cam 
z 


and  insert  into  (5,27); 


^n^T  iq  (n— 1)  a 


Si  (2/N  ) [sin  om.a  - — e ' sin  qa" 

' 2 1+g  S_ 


-n  T 


With  the  use  of  (A, 8),  we  obtain  the  desired  estimate  (5 


V is  defined 

and  obtain: 

(A.  7) 

(A.  6) 
(A.  9) 

(A.  10) 

.29)  . 
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9.  A number  of  papers  concern  the  numerical  study  of  slabs  of 

various  thicknesses.  As  an  example,  the  resolvent  method  that 
"works  for  finite  crystals  (in  practice  up  to  about  1000  atoms) 
with  two  surfaces  parallel  to  a chosen  plane":  A.  van  der  Avoird, 
et  al , Phys . Rev . BIO . 1230  (1974) . However  some  authors  have 
tackled  the  semi-infinite  solid  directly,  e.g.  J.A.  Appelbairr  and 
D.R.  Hamann,  Phys.  Rev.  Letts  31,  106  (1973)  and  3^,  .225  (1974), 

10.  See  definition  and  plots  of  l.d.o.s.  (called  LDS  therein)  in  J.R. 
Schrieffer  and  P.  Soven,  Phvsics  Todav,  April  1975  issue;  pp,  24- 
30, 

11.  See:  J.  Koutecky,  Phvs . Rev.  108 , 13  (1957)  ; J.  Kouteck;.-  and  M. 
Tomssek,  Surface  Science  3,  333  (1965);  D.  Kalkstein  and  P.  Soven, 
Surface  Science  26 , 85  (1971)  . 

12.  Well  known  discontinuities  in  the  d.o.s.  in  the  variational  and 

electronic  band  structures  of  solids;  see  van  Hove,  Phvs . P.ev . 89 , 
1189  (1953),  who  firs^ derived  them  on  topological  grounds,  or  any 
solid  state  text  such  as:  J.  Ziman,  "Principles  of  Theory  of 

Solids",  Cambridge  University  Press,  1964. 

13.  If  it  is  not,  Schmidt  orthogonalize  and  normalize  them  first. 


14.  This  method  of  solution  of  the  multiple-scattering  problem  may  find 
other  applications,  as  in  the  study  of  not-so-dilute  impurities  in 
solids.  For  another  application,  see  Section  VII  herein. 


15.  One  often  has  to  be  content  with  expressing  the  scattered  states 
as  solutions  to  a Fredholm  integral  equation,  to  be  developed  as 
an  infinite  series,  or  approximated  variationally , In  the  present 
work,  we  develop  only  the  exactly  soluble  cases. 


16.  Aside  from  the  books  cited  elsewhere,  a considerable  bibliography  is 
available  in  Physics  Today.  April  1975  issue.  The  pioneering  works 
are  by  I.  Pamm,  Phys . Zeits . Sowi . 1,  733  (1932)  and  W.  Shockley, 

Phys.  Rev.  56.  317  (1939), 

17.  But  not  all  — in  some  schemes,  the  surface  solutions  are  matched 
to  the  known  bulk  solution  at  the  pth  plane;  see  Appelbaum  and 
Hamann,  op.  clt.^ 

18.  P.  Hohenberg  and  W.  Kohn , Phys.  Rev.  B136 , 864  (1964)  and  W,  Kohn 
and  L.J.  Sham,  Phys.  Rev.  A140 . 1133  (1964) , applied  to  surfaces  by 
N.D.  Lang  and  W.  Kohn,  Phys.  Rev.  4555  (1970  and  B3,  1215  (1971), 
and  J.  Perdew  and  R.  Monnier,  Phys.  Rev.  Lett.  37 . 1286  (1976)  and 
D.C.  Langreth  and  J.  Perdew,  Solid  State  Communication  17  , 1425  (1976)  . 

19.  Discussed  in  many  texts,  such  as  Ziman,  op.  cit  (reference  12  above). 
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Figure  1 

Contours  of  constant  e for  "cubtun"  In  the  2D  B.Z.'s  appropriate  to; 

mm 

^ (001)  surface  orientation  a-(cos  *=0®  ^y®o 

( (Oil)  surface  orientation  ^-(cos  k^a^-^2]cos  ^|) 

(schmatically) . Regions  of  instability  against  surface  states  are  in- 
dicated by  shading,  and  are  obtained  from  equation  (5.11)  and  Figure  2). 
These  exist  only  for  the  TOll)  orientation  and  over  the  restricted  energy 
range  -1  < c £+1;  the  full  band  extends  frotr.  -3  to  -t- 3 . 


Figure  2 

Loci  of  constant  for  (Oil)  surface  of  cubiun:.  Shaded  region  near 

■ 0 indicates  instability  against  the  forrcation  of  bands  of  surface- 
states  for  arbitrarily  weak  surface  perturbations.  Such  regions  of  in- 
stability are  generally  curved,  for  arbitrary  band  structures. 


I 


rigure  3 

Typical  energy  contours  at  t * (n,k  ,k  ) . The  independent  variable  k 

X y z 

is  just  one  of  the  roots  of  the  equation  e^(k^)  * const.,  vhere  e,(k^) 
is  given  in  equation  (3.6).  The  total  nuaber  of  positive  real  roots 
of  this  equation  is  Q,  indicated  on  the  left-hand  colunsn.  In  the  text 
(Sections  III  and  IV)  it  is  shown  that  only  Q of  these,  coupled  with  a 
number  of  the  complex  roots,  survive  the  imposition  of  surface  boundary 
conditions.  It  is  also  proved  that  surface  amplitudes  vanish  at  points 
marked  (*)  where  Q changes  discontinuous ly . 


Figure  4 

Bound  §taces  and  their  effect  on  the  d.o.e.  (schematic). 

(a)  Graphical  construction  of  solution  of  equation  (5.6)  for  a 
single  band,  "o"  indicates  the  bound  state,  "1"  the  conjugate  resonance, 
for  the  example  of  an  attractive  perturbation  (g  < 0). 

(b)  L.d.o.s.  for  (Oil)  cublua,  m = 1 . For  week  perturbations,  as 
assumed  in  this  figure,  the  surface  (bound)  states  contribute  only  over 
a narrow  range,  indicated  in  the  shaded  regions  of  Figures  1 and  2.  In 
the  bulk  (tn  ■ •)  , the  surface  contributions  vanish,  and  van  Kove  singu- 
larities re-appear. 


Figure  5 

Principal  surface  regions  and  perturbation  (scheciatic)  . 

(I)  Vacuisn  region,  extending  from  z « —Pa  to  -L  (L^~*  •)  . Eigen- 
functions are  asumptotically  plane  vaves  having  positive  energies 
^2k2/2ffi;  a complete  set. 

(II)  Finite-sized  vacuum  interface  region,  0(Pa  ),  Two  sets  of 
functions  are  needed  for  completeness. 

(a)  Plane  waves  are  modified  here  by  orthogonalization  to 
eigenstates  of  the  solid  and  to  each  other. 

(b)  Eigenstates  of  the  solid  leak  into  this  region,  to  vary- 
ing extent. 

(III)  Finite-sized  solid  interface  region,  ©(Pa^) , also  two  sets 
of  functions  are  required. 

(a)  States  of  the  solid  are  affected  strongly  by  varying  Hg . 

(b)  Some  vacuum  states  leak  into  here. 

(IV)  Asymptotic  solid.  The  eigenstates  are  the  essentially  com- 

plete set  of  Bloch  functions;  the  solid  has  arbitrarily  large  length  Nja- 
(although  — 0).  It  should  be  nooed  that  there  la  a distinction 

between  the  vacuum  states  and  the  positive  energy  states  of  the  solid, 
regardless  whether  matrix  elements  connect  them  or  not.  Both  sets  are 
required  for  completeness.  Cross-sections  are  N^Sx^ySy , same  in  all 
regions , 
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We  ^UlJy  llie  elioets  ol'u  surface  perturlialion  on  a semi-infinile  solid.  Hie  bulk  band  structure  is  found  to  deter- 
mine wlielber  or  nol  the  surtacc  is  unstable  ajiainsl  the  torination  ol  electronic  surtace  bands.  .A  criterton.  invidvinp 
a sinele  coinpitnent  of  the  bulk  mverse-elfective-niass  tensor,  is  derived. 


I'.xperiincitUil  |1|  atij  tlteorelieal  * ' analysis  have 
by  now  eslablisltej  that  tlie  stnl'aces  of  solids  ol'len 
dil'tcr  In  eryslallograpliie  class  (.s|  as  well  as  in  electfo- 
nic  ptcipcrlics*  - from  the  bulk.  The  fact  lliat  not  all 
surfaces  have  beeit  fottttd  to  be  reconstriicted  led 
me  to  seek  a criterion,  and  ultimately  to  pose  a simpler 
problem:  uitder  what  cireunislances  will  an  infinite- 
simal perturbation  produce  a band  of  surface  states  ' 

We  shall  denote  this  an  “intrinsic  instability",  for 
while  any  arbitrary  surface  wilt  have  btwmd  states 
I'ltly  if  the  perturbing  potentials  at  the  surface  are  suf- 
ficiently great,  one  with  intrinsic  instability  it/uui  s 
must  *■’ . There  have  been  previous  atleinpis  at  sysie- 
niati/ing  criteria  for  the  e.xisteuce  of  surface  states 
|6):  the  charge  redisiribitlion  effected  by  their  exis- 
tence may  additionally  be  responsible  lot  the  observed 
surtace  reconstruction  |.7|.  In  this  paper  we  offer  nia- 
theniatical  pnuit  that  a surface  defining  the  v r plane 
is  intrinsically  unstable  iff  a component  of  ihe  bulk 
inverse-mass  tensor  a.,  vanishes  along  a set  iif  points 
or  lines  in  the  Ihillotiin  Zone.  We  thus  foeiis  on  the 
topology  of  the  hand  slrticltire  as  a root  cause  of 
the  surface  instability  : the  exact  naltire  and  strength 
ol  the  snrlace  perturbations  //,,  are  of  secondary  im- 
porlanee.  (Our  analysis  does  nol  require  any  particu- 
lar methodology  such  as  light-binding.  NTI-,  ,AI’W,  etc.) 

* I hiv  research  supported  by  a grant  I'ruiii  the  OlTiee  of 
.Naval  Research. 

Some  theoretical  concepts  are  reviewed  m |2|,  note  that 
they  number  the  first  atomic  plane  ii  = It.  w hereas  in  the 
or**sent  work  it  is  ri  = I . See  also  chapters  liy  .tones  ,ind 
■>  K m I I |. 

♦ 2 

.a  .'ll  iis  Jorinalioii  lovcis  bi>und  to  (Ik*  nei^hhor- 

IumkI  ot  a surtacc'surtaco  stales.  Seerel.  )4)  .im)  ecj,  (4)o|  (In- 
present  paper. 

4(>H 


The  prototype  material  consists  of  atomic  planes 

at  Zy  = m^a.  with  iiij  = 1.2 V_.  the  plane  index.  The 

interplanar  distance  a is  a tunction  ol  the  orientation 
of  Ihe  surface  relative  to  the  crystal  axes,  of  course. 
The  unperturbed  Hamiltonian  //|,  olTlie  lerniinaled 
solid  has  eigenfunctions  (r).  in  which  k,.  and  A,, 
are  the  usual  cry  stal  momenta  but  A,  is  restricted  to 
the  interval  0 < A-<  n'u.  I hese  eigenfunctioiis  have 
been  shown  |7|  to  lake  the  Ibriti: 


' . ' ' * 
in  vvliicli  R^)  is  ilio  usual  Wannicr  riinction 

(the  l■o^tricr  transi'orm  of  the  Ifloch  functions  foi  the 

nth  band  using  periodic  boundary  conditions).  While 

terminalion  ol  the  solid  preserves  the  translalional 

invariance  in  the  .v  c plane  (witness  the  exjioiienlial 

factors)  it  has  the  effect  of  replacing  .\2  ' “expliA-/  ) 

by  the  components  c>t  a new  A2  X .\ . unilaiy  mat t i\ 

F„k(Z/).  I he  construclioii  ol  this  ftinction  and  its 

analytical  properlies  ate  elahoitiicti  in  a conipanion 

paper  |7|.  We  iptote  one  lesiih  suj-i|iose  llie  eigenvalue 

crtnesponding  to  ( 1 |,  denoled  c„(A^ , A,. . A . ) has  a 

mininiunt  or  ntaxinumi  at  A.  = A,„  (for  fixed  A^..  A, 

aiul //).  Defining  (S  (A. 

iniporlance. 


A„j  )'i.  the  lolkwcing  ts  ol 


li"'IF„^(w,i/)|- 

h *{) 


i (2  .\ . )i>rh ■ 


(21 


lor  ut,  any  (loMlive.  finile.  Inlegei . liy  cotiliasl.  deep 
in  the  bulk  (ni^.  - ■f'»)  T ,|U“  ~ ^ - is  a constant  in. 

tiependeni  ol  A . and  ouiside  the  solid  (m,  - (I  oi  ne- 

l trie  1C  tv c. luce  cull. ICC  pci  lurlvuiuitc  muci  .ilw.o  c cxici  I ,ii 
cs.impic.  the  M.idcliiiii'  poiciiti.il  miici  dcu.iic  m Ibc  iicicb- 
bothood  ol  .1  ciirl.uc  Irom  pc  bulk  limitmc  v.iluc.  See  cucb 
klt'sk  usMxwiN  ,>s  >n  I 5 I 
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ll\o  t I ,,1^  D-  I J ) IN  pilicK  .1  proiX'i  1>  ol' 

iIk' mtciiiii  ik’I;JiI'ui1uhh1  ul  llio  sin  1.k\'. 

Wo  nou  inliixluoo  a poi  liii  Iniii:  I l.miilloni.iii  //^ 
uIikIi.  loi  MiiiplioiK  . IN  taken  diaiinnal  in  the  Wanniei 
Innetiim  K'l'iCNcntation  '■* . It  InllowN  tlial  its  main\ 
Ntriiv  Uno  in  ilic  baNiN  set  ( I I is. 

n'A  .1 


S .A’,,.,,  ( I /'.t » f ( /'cl  A . A l.'S,,.,,  ■ . I ■' ) 

/’  ■ 

li\  Liinstnielnin.  //^  jirosoives  liaiislational  invaiianee 
in  the  V r piano  aiul  thus  oonsoi\os  numionia  A\..  A,  . 

I ho  opiaiuitios  y,,  tlio  ooiipling-oc,mstanls.  paiainoln/o 
llio  sironatli  ot’tho  portin  halion  at  liio /Mh  piano  ror 
tlio  //  th  band. 

Wo  now  stnd\  bound  states  ol  //,,a  //^.  ibo\  are  ro- 
i|iniod  to  bo  of  the  I'oiin; 


VAv.A,..s"\-  ' ‘ p ^'-i.Av.A, •'*'hA*''*  ■ 


(4) 


Ibis  Is  to  bo  oontiaslod  with  bulk  stales,  wbiob  lake 
ibo  fonn: 


VA,.A,.A;  ‘ p''’Av.A;v'''A,a, 


X/„A«A  '‘'»A 

III  the  piosonoo  of  the  sint'aoo  |>oiliirbation  //^.  I bo 
oxaoi  soattomii;  iboorotio  oalonlalioii  ol7,,,^.  is  givon 
m a sopaiato  paper  . I\n  present  purposes  it  is  <.piite 
siilfioionl  to  siinl>  the  bound  stales  (4).  if  an\ . aiul 
ibeir  energies  /•./(, a^.. A,  .s-  -Ws  A^..  A,  and  n are  kept  fi.xod. 
we  omit  tboso  snbseripis  in  wbat  follows  lor  greater 
topograpbieal  elaritv . I list.  SebrfKliiiger's  eijiiation 
(//ii+  \ ields  the  folbtwing  Oiiiiaiion  for 

the  eoelTieients  /-'(A. ): 

Zy/•■(A.)t(A.)'l'^7  Z/  V 

A_-  ' ■ - A.-  - ■ ■ 


- / \ p / (A,  I'l7_.  . 


((') 


in  wbieb  e ( A . I is  the  bulk  band  siineiuio  t„(  A I w itb 
stipei  llnoiis  iiidiees  siippiesserl.  I ipialiiig  eoei  lieients 
of  'b(i,'s  we  leadiK  obtain  a set  ol  liiieai.  boinogenoous 
equations  loi  the  /■  (A,  fs.  and  a eondition  loi  a non- 
ti i\  lal  solution  to  exist 


De'tl'A,,,,.  J 0.  (7 

1 be  anaK sis ot  tins  seeulai  deteinnnanl  is  the  main  oh- 
|Oet  ol  tins  papoi . 1 be  .V^,  aie  intogials 


) = ■'  fdA. 


«A</4U/,  _.(/'■  1 
/ t ( A ■ ( 


We  bare  absoi  bod  a lae  toi  of  ,\  ! " in  / ^ to  loi  in  the 
ti/^  independent  ol  \’_  in  the  tbeiniode  nainie 

(bilge  .\  I limit.  I or  / ^ wiibin  the  eontinniim  lange  ol 
( (A . ) the  .S'^,  aie  eoni|'lex  and  tboieloro  ( 7 1 has  no 
sobilion.  I bus  tile  bound  stales  be  oiitsule  the  extrema 
at  A,„  . w bile  cpnie  possible  oeei  lapping  some  ol  liie 
bulk  eon  t inn  mil  of  states  foi  dil  leient  (A  ^ . Aq  )'s  or 
iliffeienl  band  index  n 

\s  we  have  seen  in  ( 2 ).  the  nmneialors  in  ( H l vaiiisb 
when  A,  approaebes  an  exiieninm.  riieiefoie  tlie.S'^, 
are  ortlinariK  finite  mtegials.  no  mailer  bow  elosek 
/•^approaebes  the  eontimnmi  edge  at  e(A„,  I.  Now  let 
us  examine  i (A. ).  i.e.  t„(A  ).  iieai  an  extiemum. 

t ( A I “ e^jl  Aq. . A . Aq^j  ) + t rV-_(  Aq. . A . Aq^^  It  A _ Aq^j  l"i/~ 

+ ()((Aq_  A,„|4),  ('ll 

where  a,_  “ d-rq,(A  I il(uA- 1-  is  the  rr-eomponenl  of 
the  inverse-efleetive-nuiss  tensoi . to  be  eeabialed  at  the 
point  (A^..  A^,.  A„,  I in  the  utli  band  * * . W hile  noinialK 
positive  at  a niiiiinunn  and  negative  at  a maximimi.  this 
t|uantit\  eoiibl  be  /eio  at  eilbei , eaiising  .S'^,  1"  rb'eige 

as  l/-.'^  i4A„|l|  ' ■' w ben /' ^ appioaebes  ( ( A I.  1 bus. 

the  erilerion  I'oi  ( "I  to  have  a solution,  when  the  eoiipl- 
iiig  eonstanisg^,  aie  .iibitrarilv  small,  is  ineielv . 

a._-(Aq..A,..Aq„l  = ()  . (lOl 

W'berever  ibis  is  salisfierl.  we  neeessaiilv  have  snrfaee 
states.  Ibis  nia\  oeeiir  ai  diseiele  poinis  oi  along  eiiives 
III  i be  (A ^. . A,  I plane.  Ik  eonimmix . ,it  finite  values 


In  a mure  oniipleie  stmb  nl  ilie  I’erliul'ed  sutl.ue.  iiinlei 
rpep.it.ilion.  we  consKler  inler-txuul  iiulns  elements,  solve 
l.ir  se.illenm-  states  and  iliaree  redisinbnn.ni.  densilv  -n|- 
slates  luiu  lii'iis.  elv. 


I he  Iveh.ivjnr  ol  ( , eq  I 21.  is  puii-lv  eeoineliKal  .nut  nn- 
rel.Ueil  to  ihe  in.ienitnite  ol  u - - I Inis  .niv  diveieeni  beli.ivnn 
ot  the  .S'^,  inlevi.its  IS  e.nised  solelv  1>\  the  eneiev  denonii- 
n.ilors. 
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I'l  a liniie  neigliborluuKl  of  tliese  poinis  or  curves 
also  rtoccssarily  sauslics  lire  condiiioii  (7)  for  llic  cxis- 
iCMce  ol  a bouikl  surface  state. 

Us  way  ol  an  illuminatinit  application,  consider  the 
s.c.  tight-binding  batid  structure:  e(A)  = A'(cosA:^.«  + 

cosk^M  + arsA.i/l.  for  a phy  sical  surface  perpendicular 
to  the  ( lOU)  axis.  The  band  extrema  are  at  = 0 or 
IT  u.  and  takes  on  the  values  +K  at  the  one  and 
K at  the  other,  lor  all  (A^..  A^  ).  Thus,  the  crilerioit 
0(.|.  I 10)  for  intrinsic  instability  is  met  nowhere  in  the 
do  Urjllouin  /.one.  On  the  other  hand,  tilting  the  sur- 
lace  perpendicular  to  the  (Oil)  direction  results  in  new 
lattice  spacings:  u,.  = u,  = u 2*  - and  = a.  and  in  the 
new  coordinate  sy  stem  (A_  axis  must  always  be  perpen- 
dicular to  the  physical  srrrl'ace)  the  etrergy-batid  strirc- 
titre  is:  e(A)=  A(cosA^.«  -1-2  cos  Aj.r/^,  cos  A-u, ).  The 
extrema  are  at  A,„a  = 0 or  ,t,  and  eq.  ( 10)  now  has  turn- 
iribial  solutions  along  the  segments  k^.Uy  = rr,  for 
all  A^..  This  orientation  surface  is  intrinsically  unstable, 
whereas  the  other  orientatiori  was  not,  for  the  saute 
material. 
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Klectron  states  in  random  alloys  with  short-ranj»e  order 

Paul  lOoom*  and  Daniel  Mams’ 

Heifer  (iruiiuuw  School  nj  Sacncf.  )\‘\hi\u  rMiwr.jfy.  \V« 

(KccciU'il  < Auj!U'.( 

Ue  present  an  aeeurale  and  eeommneal  ileralive  methiKl  nf  ealelilalinp  the  ener^N  lerels  ol  i disordered  or 
parlls  ordered  random  alloy  Results  presenteil  tor  one-  atitl  ihree-ditnenstonal  snnple  euHie  latti.  es  tonipiire 
I'asorahK  with  e.xaet  ealeulatiotls  We  also  present  the  systetnatie  elVeets  ol  partial  short  range  order  in  three 
dimetistons  .-\  theory  of  the  one-parttele  propagatots  is  presenletl.  and  the  theory  of  eleelrnal  eomlitt. tiuty  is 
deveK'ped  in  the  eontext  of  iiur  new  metluHl  (^tir  foltiuilas  satisfy  the  exaet  eonsetxation  laws 


I IM  RODl  ( IION 

The  study  of  electronic  and  vibrational  spectra 
of  disordered  alloys  is  currently  one  of  the  prin- 
cipal concerns  of  solid-state  physics,’  stimulated 
by  the  nutstandiny;  successes  of  the  colie rent- 
potential  approximation"’^  (CP.A).  now  ending:  its 
first  decade.  Although  efforts  to  improve  our 
understanding  beyond  the  CPA  have  not  all  met 
with  the  same  good  fortune,  there  have  been  re- 
cent exceptions.  Cluster  methods''  ’ have  been 
devised  which  are  accurate  enough  to  reproduce 
the  "peaky"  structure  of  the  density  of  states  p(x>), 
whicli  they  sometimes  do  (notably  in  one  dimen- 
sion') with  startling  fidelity.  W e have  been  work- 
ing along  such  a cluster- type  approach,  and  have 
found  an  extremely  simple  method  translaling  di- 
rectly into  a computer  algorithm.  While  unsuited 
to  the  theoretical  study  of  Lifshitz''  tails,  our 
method  has  permitted  us  to  reproduce  many  of  the 
other  known  results  over  the  theoretically  per- 
mitted range  of  energy.'  even  near  tlie  energy  max- 
ima and  minima,  and  additionally,  permitted  to 
study  of  the  effect  of  sliort- range  order.  Along  with 
Lifshitz,  we  envt.sage  tails  in  p(x’)  at  the  energy 
maxima  and  minima  as  arising  from  accidental 
correlations  in  increasingly  largo  clusters,  of  a 
size  that  for  practical  reasons  we  are  not  at  pres- 
ent capable  of  handling;  however,  the  simplicity 
of  the  present  method  may  suggest  a natural  ex- 
tension to  cover  this." 

I'he  basic  outline  of  our  paper  is  as  follows: 

In  Sec.  II  we  present  a method  for  the  calculation 
of  the  single- l)ody  Circen's  function  in  the  presence 
of  an  arbilrarv  number  of  inipurilies.  We  tlien  dis 
cus.s  liow’  our  procedure  can  bo  implemented  by  tlic 
use  of  a convergence  factor  il.  Section  III  is  de- 
voted to  an  analysis  of  Ihe  meaning  and  uses  ol  lh< 
complex  sell-energy  'i  witliin  llie  conlexi  ol  .1  dis- 
ordered medium  llesulls  1 rom  our  method  .ire 
presented  in  Sec  IV.  Including  lhi‘  cticcis  ol  shorl- 
r.mge  order.  Hevond  this  m Sec.  V we  m.ike 
further  approxini.ltions  that  .illow  us  to  deter- 


mine Lf  i(s’).  Section  VI  is  concerned  with  the 
development  ol  a transport  theory  compatible  with 
(T,  along  the  lines  of  Baym  and  Kadanoff.  ' 

II  ( LI  SI  I R (,RI  I N S IT;N(  IION 

Let  the  Hamiltonian  for  the  electrons  within  a 
single  light-liinding  band  in  a hypercubic  lattice  in 
1)  dimen.sions  be 

ui  1 

with  Tj^r.{2l)Y'  for  /./  nearest  neighbors  and  zero 
otherwise,  ii)  the  W mnier  state  at  the  lattice 
point  and  I’,-  the  p 'ential  which  takes  on  one 
of  two  values  dcpendii  hether  atom  A or  B oc- 
cupies the  i site.  We  construct  the  resolvent  oper- 
.ilor  t.(z)  and  its  various  matrix  elements: 

(;(.-)  x|  z (i’_ , (2) 

in  which  we  reference  tlie  nper.itors  to  ,1  complex 
"optical  potential"  !:!(-')  merelv  ,is  .1  device  to 
enhance  the  convergence  ol  subsequent  exiian- 
sions,  with  : equal  to  the  Ireqiiency  w.  extended 
to  the  complex  plane. 

For  lliose  readers  lamili.ir  with  the  Cl’.-V,  it 
is  important  to  note  that  our  new  dcp.irture  con- 
sists principally  in  dissociating  the  complex  self- 
energy parameter  I'lxf  from  Ihe  M(e-d;.ig«nal 
averaged  Gri'en's  function  W hcreas  in 

Cr>.-\,  knowledge  of  tlie  one  implies  the  other,  via 
tile  relationship 

1'-! 

('  I’  A 

our  experience  indicates  th.it  it  is  lietter  to  tre.it 
ll(x')  merelv  .is  a convergence  |ia r.imeter , one  to 
he  chosen  as  ,in  111/  lin,-  ,iid  m Ihe  c.ilciilalions 
r.ilhor  Ilian  bv  tedious  .ind  iiniiecess  iry  self-coii- 
sislenev  coinliltnns  As  bv  Kq.  (2)  the  ('luc/  (■„„(-') 
.ire  .ill  unit  |H'iidenl  ol  I.  in  .iny  incxni/r  np- 
'i.iKiii  n'liih  to  I we  have  l.ititiide  in  our  choice 

ol  ",  i i.  ,1-  iliscussed  below  ,iud  w e pick  I he  sim 
oil  St  po-...ibli  ■ I I lo|-  which  our  c.ilcuhllcd  (.  is 
i)i)ii  .xiniiiclx  sl.ilioii.irv. 
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We  lU'xl  deliiie  ,i  modified  resolvent  openilnr 
npiiropruite  to  ihe  cnse  m whieh  one  sets  1', 
^0.  wliore  we  define  f,  '(  ,i)  nnd.  in- 

die.iting  tile  elimination  of  the  loealized  lluetiia- 
tion  potential  at  this  site  liy  ( we  have 

The  full  resolvent  (2)  ran  be  expressed  in  terms  of 
the  modifieation  in  (-1)  by  the  use  of  tile  operator 
identity  (.1-  «)-‘  - .-V' + .-l-'/Jf  ; 

(■(^)  (;"V)  + f;‘"(-’)(r-  -;),(.(-)•  (5) 

Because  the  perturbation  is  diagonal  intheWannier 
representation,  the  matrix  elements  are  easily 
found: 

(G) 

For  the  calculation  of  the  density-of-states  func- 
tion/i(w)  (-1  it)  liii(;„.,(j.' i-/e)  only  the  configura- 
tionally averaged  Green’s  function  Tit)  is  re- 
quired. For  the  one-particle  propagators  Gjj 
the  .iveraged  Fourier  transforms  of  .ill  + 

are  needed.  Kquation  (6)  is  now  iter.ited.  De- 
fine to  be  the  modified  resolvent  oper- 

.itors  with  the  fluctuation  potentials  f at  sites 
/ .ind  I removed.  By  a repetition  of  the  above, 
we  have 


,’U)  = (■ 


( I »J  > I • / ' ( I , ■ J 


n j ' ' 

(•7) 


The  in.itnx  elements  decav  exponentially  w itli 
dist.ince  thus  the  e.xpansionsfG)  and  (7)are  in 
.1  symbolic  "parameter"  i defined  as  (.‘„y  l'^. 
which  IS  "small"  for  small  f,  and  "e.xponentially 
small"  at  l.irge  I'j.  I'lu-  /irocevM'.s  i'>i  <iinl  i,  i (iiv 
lt>  /)(’  mi\  niiiiihi’r  oi  //)//(•.^.  mi/il  llir 

l>y(uliciil  rliistry  >,tir  /.•>  ncli/i’i'i’i/ J Ter- 
mination. by  truncation,  ol  the  series  consists  of 
approximating  the  most  distant  (.'s.  i.e..  those 
with  the  largest  number  ol  .superscri|ils.  by  their 
value  in  the  average  oplical  potential.  Thus,  il  we 
stop  at  (7).  the  approximation  consists  in  re- 
pl.icing  (.'.I"  ^’(.')  by  e 1 1 - (7"  i il)  |"'  ni.i  The  con- 
figurational averages  over  .ill  the  explicitly  re- 
tained fj  are  then  performed,  and  .ill  G's  obt. lined. 


Ill  t IKIK  I Ol  1 

We  now  come  to  our  princip.il  point  ol  departure 
from  other  methods — our  choice  of  Our  results 
would  depend  cruci.illy  upon  ')  except  for  the  fol- 
lowing observations.  Since  Ihe  beh.ivior  of  the 
local  cluster  is  the  dominant  char.icteri.stir  of 
disordered  svsiems.  we  expect  results  insensitive 
to  Ihe  particul.ir  choice  ol  if  the  cluster  size 
IS  sufficieiil Iv  l.irge. 


We  require  .1  simple  fuiiclion.il  lorm  lor  . tli.it 
.illows  lor  states  out  to  the  bands  limits.  This  ex- 
cludes the  use  ol  which  is  known  to  jiroduce 

bands  th.it  are  .ilwats  too  narrow.  We  restrict 
the  range  ol  possible  b's  by  re(|Uiring  that  it  obey 
dispersion  relations,  insuring  that  our  approxi- 
mate (,  is  analytic.  Furthermore,  a functional 
form  is  desired  in  which  (.  is  accurate  in  both  the 
weak  as  well  as  strong  scattering  regimes.  Be- 
cause of  the  local  nature  of  highly  disordered  sys- 
tems. our  choice  becomes  more  critical  tor  small 
potential  differences  where  elfecls  are  more  ex- 
tended. Our  input  is  the  iiiiD  which  we  lake  as  one 
or  more  step  functions,  nonzero  only  within  the 
theoretical  band  limits.  Bei;  is  then  determined 
from  the  following  dispersion  relation: 


s./'. 


*-/ o 


(8) 


This  is  sufficient  to  make  our  .ipproximation  to 
(-'(z)  satisfy  causality.  The  density-ol-stales  sum 
rule.  ( (fw  (ifw)  1 IS  itself  a beneficial  con- 
sequence of  the  aiialyticity  of  our  approximate 
G(')  and  its  resulting  1 z dependence  in  the  asymp- 
totic limit  as  we  discuss  elsewhere."  We  verified 
that  in  all  cases  studied,  the  sum  rule  on  piv) 
was  satisfied  nuniericallv.  If  we  choose  the  con- 
stant Ini')  to  be  ol  magnitude  of  IniD^j.^.  then  in 
Ihe  weak-scallering  and  low-concent  r.ition  re- 
gimes (,{^)  will  be  quite  similar  to  so 

that  .iccurate  results  ran  be  e.xiiected  in  .ill  re- 
gimes. 

Before  we  proceed,  the  way  in  which  we  use 
IiiD:,  must  be  more  clearly  outlined.  In  Ihe 
.iccompanying  Fig.  1 we  display  the  two  basic  be- 
havior p.itterns  of  Ini'll, as  observed  by 
Velicky.  Kirk)iatrick.  and  Khreiireich.'  It  should 
be  noted  that,  here  also.  RelJ  and  liii"')  are  re- 
lated by  the  Kq.  (8)  has  to  describe  every- 
thing in  the  ri’,\;  it  determines  b.md  gaps,  pe.iks 
in  the  density  ol  slates,  .md  the  gener.il  overall 
scale.  Most  of  these  results  (e.g..  band  gaps  and 
complic.ited  structure)  are  better  obtained  by  our 
detailed  calculations  of  the  correlated  scattering. 
We  hypothesize  that  the  most  useful  information 
from  Cl’.\  is  conl:iined  in  Ihe  general  overall  mag- 
nitude of  Imilj.,.,,  Oiierationally.  in  Fig  1(a).  we 
would  ignore  values  of  Im';^.,,^  ^ from  the  region 
of  its  maximum  as  well  as  the  extrenieties  of  the 
band  In  the  former  range  of  energies,  we  ex- 
pect ('xceplional  scattering  because  it  is  easiest 
for  these  stales  to  make  l r.insil ions  due  to  band 
overlap,  whereas  at  the  band  edges  the  spectrum 
will  be  least  disturbed,  according  to  the  s.iiiie  con- 
siderations .Viiv  value  from  Ihe  shaded  region  is 
then  acce)ilable.  In  terms  of  parficle  lifetimes,  we 
will  obt.iiii  Ihe  l.irge  .ind  small  t r.insil  ion  r.ites  be- 
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KRi.  1.  Ufal  atui  iniai;inar\  paits  t>[  the  eomplex  seli- 
■ni  iLtt  in  the  Cl’A.  rhe  Iwo  sets  oteiirves  are  indieative 
ot  till'  tV|X'  of  results  that  ean  he  expeeleil  from  this  ap- 
l-iroximalion.  (a)  A situation  where  th('  alloy  haniis  over- 
lap; (hi  the  ease  of  split  hands.  Ktierev  i.s  in  units  of 
half-handw  idth.  This  is  an  adaption  of  a fitfure  from 
Kef.  1’. 


cause  we  almost  solve  Uie  eiftenvalue  problem  ex- 
actly for  each  configuration  and  this  is  clearly 
e(|uivalent  to  a perturbation  aiiproach.  As  for 
concentration  dependence  in  ^cpa-  ''ill  obtain 
correct  behavior  .sim)ily  because  we  wetph  each 
confisjuration  by  its  appropri.ite  probabilitv.  Thus 
we  are  able  to  include  both  the  dynamical  and 
statistical  .ispects  of  the  problem. 

In  Km.  1(1)1.  the  s.ime  an.ilysis  leads  us  to  ipiiore 
the  very  lari;e  values  ol  Im  1(  pa,  both  subbaiids 
Here  tliout;h.  the  maitnitudes  are  consideralily  dil 
fereni  le.idinit  us  to  suspect  that  two  diflerent  con- 
st.ints  .ire  needl'd.  Kiirther  det.iils  of  this  case 
will  be  eliicid.iled  in  the  followini;  e.x.imples 


IV  \\  UA  SIS  Of  Kl  St  1 I S 

We  lirst  consider  the  canonical  one- dimensional 
tipht-bindinp  binary  alloy  for  three  different  scal- 
terintt  strengths  at  a 50-50  concentration.  Figure 
2 compares  the  results  of  one-,  three-,  and  five- 
cluster  calculations  for  p(w)  when  K,  0.5  with 
exact  results.  We  see  in  this  example  the  develop- 
ment of  tlie  peaky  structure  associated  with  special 
clusters  of  atoms  as  our  cluster  size  increases. 
Proceeding  to  a larger  scattering  strength 
(r,  = il.O).  we  expect  that  the  local  configurations 
will  play  a more  prominent  role  becaii.se  of  in- 
creased wave  function  localization.  As  shown  in 
Fig.  3 we  successfully  reproduce  most  of  the 
structural  details  of  fifwl  for  a five  cluster.  To 
check  the  degree  of  insensitivity  in  our  five- cluster 
model  we  varied  llmll;  within  the  limits  given  by 
IniDcpA  AKd  found  little  cnange  in  the  overall  pat- 
tern as  shown  in  Fig.  4.  This  indicates,  numeri- 
cally that  the  resulting  C is  stationary  and  that  H 
is  optimum. 

The  scattering  strengths  are  now  increased  to 
lb  - -•  2.0,  providing  a critical  evaluation  of  the 
methods  capabilities  (larger  scattering  strengths 


l HI,  g.  t'omp.-iri -on  ol  ono-  , Ilirn'-  , ;ind  liw  * Im-li  i 
calc'ilation.''  for  p(  i ii.xing  Imyy,.,^,  u.l  >,  with  f, 
ii.V,  I-  11,  1.  Itnokavouipi  (hixtogr.-iml  i:-  ox.n  1 i*  -ib 
o|  Kol.  I.  I.iiorgv  is  ill  iiiiils  of  h.'ill  haiidw  idlh. 
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l‘Ui.  .■>.  Cotiiparison  (>1  otu'-,  tliroc-,  rind  n\'('-clustc‘r 
nalfiilatiims  lor  p( , ) in  on«’  ditnonsion  usinr; 

= n,r>,  n ilh  l',=  1.0andr  (),.*).  llark|rround  (hisloi;rain) 
is  oxaot  ri'sulls  ol'  Urf.  J.  Kmwxi  in  units  o!'  half- 
liamKvidth. 


are  m a sense  too  easy  because  uave  function  lo- 
calization makes  a cluster  calculation  more  plaus- 
ible). Usinn  the  exact  scattering  off  all  configura- 
tions of  five  atoms,  the  highly  discrete  spectrum 
is  well  reproduced,  as  seen  in  Fig.  5(a).  In- 
creasing the  cluster  size  to  seven  atoms,  keeping 
the  convergence  factor  the  same  as  m Fig.  5(a). 
improves  our  agreement  with  the  exact  results  as 
sliown  in  Fig.  5(b).  In  Fig.  (!  we  display  the  re- 
sults of  again  v.irying  [Im'ij  within  the  limits  dic- 
tated by  Ini'J^.,,,;  the  major  details  .ire  .igain  seen 
to  remain  stationarv.  We  have  lound  empiric.tllv 
that  'f  Imi)  I IS  too  sm.ill.  the  resuii.mt  densitv  ol 
states  IS  too  "peakv”  .md  .is  such,  repri'sent.iiive 
ol  .1  molecular  cluster,  mste.id  ot  the  solid  st.ite, 

11  lm')|  is  loo  large,  then  the  cenlr.tl  sue  pre- 

dominates. as  is  correct  only  in  the  extreme 
".itomic"  limit  when  potenti.il  I luctu.it ions  greatlv 
exceed  the  bandwidth  One  c.m  see  this  Irom  I ig 
tl  since  the  sh.irper  cirve  is  associated  with  the 
lowest  value  of  Im"  .ind  vice  versa 
In  three  dimensions  the  obvious  cluster  size  is 
seven  sites.  Figure  7 comp.it  es  our  calcul.it  ion 
with  the  Monte-Carlo-tvpe  numerical  resiilis  ol 
Albeit  /•/  .\  constant  lm\  gave  poor  results 

in  tins  c.ise.  but  the  CI’.V  c.ilctil.it ions  imniedi.ilelv 
showed  us  whv:  Im'  , w.e.  more  Ih.in  me  .rder 


of  magnitude  smaller  m the  m.ijoritv  subband  th.in 
m the  minority  -ubb.ind.  Conseiiuent  1>  we  changed 
Ini'i  to  the  step  lunction  shown  m the  ligure  lan- 
ing  the  parameters  (magnitudes  ot  the  steps)  ag.iin 
guided  bv  Cl’.A.  The  results  now  agreed  well  with 
the  ex.ict  coniiiutalions  and  wen-  insensitive  to  the 
precise  value  ol  our  parameters  as  is  evidenced 
by  Fig.  8 in  which  .i  three-step  function  was  used. 

To  illustrate  entirely  new  applications,  consider 
effects  ol  short -r.mge  ordei'  on  this  s.ime  .illoy. 
With  o the  Cowley  short- r.mge  order  parameter — 
ly,  and  1 - the  rel.itive  concentrations  and 
^‘ak  probability  of  Imdmg  .ilom  .1  .tt  .i  given 
site  when  a II  atom  occupies  a specifii'd  neighbor- 
ing site— we  have  r c^o.  /',.^^c„(l  - .»). 

' .,(1-0).  and  ‘ h''  a^-  In  llel.  13. 

o - 0,  Fore,  0,1.  (>  can  v.iry  Irom  -Oil  to 
• 10;  neg.ilive  a is  associ.iU‘d  witli  enlianced 
tendency  ol  .-I  atoms  to  be  surrounded  by  /;'s  (i.e.. 
"antiferromagnetisin").  positive  iv  indicates  en- 
hancement in  the  prob.ibilitv  of  either  species  being 
surrounded  by  atoms  ol  its  own  kind  (i.e..  "ferro- 
magnetism"). L'smg  the  same  convergence  pa- 
rau'eters  as  m our  calculation  at  o = 0.  in  F ig.  9 
we  find  distinctive  features  in  the  minority  subband 
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\|1,U  l"l'  a live  elusler  caleMlal  leii  . i|  , .ine  .ilin.-HsielKtl 

allev  w nil  f,  = 1. II  all'll-  II...  rill-  hai|ti--l|M-aksare 
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nil.Mtinn  to  (i:ilc.  lint.  I.  l-:tu-rn\  is  in  units  nt  halt-luin.iw  i.ll't. 


dt'n.sity  of  sl.itf.s  tlial  we  interpret  in  terni.s  ol  nii- 
nority-atoin  chisterim’,:  the  .sinyle  peak  ol  .> 

- 0.07  reitister.s  the  iinlikeliliood  ol  I'iiulini;  two 
.\  aloins  as  nearest  neiplibnrs.  and  Ihe  double 
peaks  of  iv  0.7  represent  the  tendeney  ol  the  s.iiiie 
.itonis  to  form  pat rs . tripU'ts.  etc.  However,  due 
to  the  sparseness  of  .\  atoms,  triplets  ,md  huflier- 
oider  ekisters  ,ire  stattstiealk  insu;ntlic.m(  lor 
these  V. tines  ol  iv. 

V lilt  IKON  I’KOI’V.MION 

So  l.ir  we  h.tve  vieveloped  a method  lor  t .ileiil.itnu. 
Ihe  site- diapon.il  eonfip.nr.ition  .iver.iyed  lireen's 
fnnetion.  We  h.tve  not  mdie.ited.  how  we  would  e.il 
fill. lie  tin  noti-sile-di.iponal  prop.u..ilors.  One.d- 
lern.ittw'  IS  to  deudop  .i  etuslei-  imdliod  lor  'he 


bitter,  similar  to  Ihe  melh.  d we  u.sed  lor  tin  lor 
liter,  .\nother.  simpler  t .iiyh  less  ,ie,  ur.iit  , .1 
lern.llive  will  be  emph.'.ed  We  lirst  ,li  Inn  ,i  m 
sell- enery.\  \’(  ' b\  the  eipl.ltion 


1 

\ • ' 


1 


rb-  rel.ilie  t.ip  . . luimerif.ilK  in\i  fled  to  oblai 
\‘  IS  .1  (tim  lion  .'I  , yji-t  or  iium.  ne.ilh  e.ib". 

bit  ed  I . J J.  I r - .11  the.  lie  III  It  Ml  >t  \ * ' .■  t w e obi  a 1 
\ • I I in  I hi.  eitt  in  I '".ple\  pi. me  .is 
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l Ui.  !».  U'lisilv  nf  slatt's  for  V^-=o.7r»,  = in  a 

rh{Vt'-<lim«-nsion:il  sinipJr  inibic  Iritlifo  with  coiua'nlralion 
(Nm  U*v  short-rango  order  paranuMrr  (I  is 
n.-b  and  0.7,  ri*S|K*ctivi*lv.  Iml'  is  lht‘ sanu*  as 
in  Fi^.  7 (dot-dash  liiU')  am!  Uel'  is  obtained  theivl'rom 
bv  use  of  K(i.  ('<!.  KncMi^v  is  in  units  half-!)nnd\Mdth. 


hence 


^ t; ' 


(12) 


as  the  definition  of  the  off- diagonal  4>lP)m>nts.  Tlu‘ 
propusiators  docay  rapidly  with  distanre  ft,„  so 
both  the  one-,  two-,  or  three-point  curve- littini; 
procedures  will  proh.ibly  sive  reasonably  equiva- 
lent results,  ,Now,  all  the  information  contained 
in  our  previous  numerical  work  is  stored  in 
the  complex  proper  self-enerj^y  part.  It  is  of  in- 
terest to  compare  1mA*  with  Im  in  order  to 
see  how  they  differ,  riiis  is  done  in  Kist,  10  for 
the  one-dimensional  alloy  of  Kins,  5 and  tl. 

In  summary  we  have  presented  a relativelv  sim- 
ple method  for  cali'ulatiiuc  the  eini'iiv.ilue  .'qiecirum 
of  a disordered  system,  one  that  avoids  all  the 
computational  pitfalls  of  self-consistent  methods. 
This  quasi- invariant  theory  is  not  onlv  hinltly  ac- 
curate, but  also  allows  tlie  liounds  on  the  frequency 
spectrum  to  be  natur.illy  determined  bv  the  corre- 
lated scatterimt  of  a loc,il  nroup.  We  now  discuss 
transport  and  develop  ,i  formalism  that  .illows  our 
numeric, il  output  to  be  used  m apiiroximat  i > is  that 
conserve  particle  number  ,ind  enemy 


VI  tl<  WSI’OKI  IN  IIISOKIII  Kl  II  SV  M I MS 

The  linear  response  of  the  current  to  the  electric 
lield  defines  lh<‘  conduct ivitv , which  we  I, ike  to  lie 
the  same  ,iloni;  the  three  principal  directions  in 
our  siniph’  cubic  structure,  Kollowmit  Velickv'' 
we  fuve  111  our  ,sinnle-p,irticle  model 


‘ f '''  (”  - '/)/’■*  • 

(13) 

where  e,  the  electric  charge,  is  unity  and  ' is  the 
fermi  function.  The  liracketed  term  is  short  hand 
lor 

({f,(.X-//)/qfi(\ -//)/.,> 

^ «olh(x  - //)/>, 6(X  - //)/>, fa)).  (14) 

a 

/),  is  the  momentum  operator  aloiip  an  arbitrarily 
chosen  principal  axis  and  ( ) denotes  confittura- 
tion  averattiii';,  Kxamin.ilioii  of  K'q,  (14)  reveals 
that  we  require  the  two-p,irticle  correlation  func- 
tion 

1 • « M (•'  1 - //)■'  UX"'  I - //)■'  • 

(15) 

Wo  can  rolalo  (G ')  lo  C by  Uu*  equation 


S 


I {(J.  b*.  liijfi  is  Im  loj  Ihi-  !)»)«•  ibnn-Ji- 

I'MKil  n|  ] . .*  iinil  U in  tin-  li\r  rliisU*t 

mnt  mu  while  1 h*’  lin*-  i.-.  I he  eni  j . - -.poniliti^ 

Im'  )l.  rniTL;\  is  in  nnil-  i*l  ijnlf  lem^lw  idlh. 


Ui  to 


I’  \ I 
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which  defines  the  index  structure  ol  the  vertex 
function.  This  equation  is  exact  for  tlie  exact  (1 
and  we  will  use  it  to  define  when  «e  have  an 
approximate  sins;le-particle  (ireen's  lunction.  \'.e 
can  place  restrictions  on  possible  verti'x  function- 
by  requiring:  the  conservation  ol  charite  and  ener- 
gy in  the  presence  of  a long- wavelength  distur- 
bance. This  leads  to  the  introduction  of  a new 
operator 

= //)•';.  tlT) 

It  is  easy  to  show  that  A must  satisfy  the  (ollowinu 
Ward-type  identity: 


which  implies  a connection  between  the  one-body 
operator  (',  and  the  two-body  Green’s  tunction 
.\lso,  one  can  show  with  the  above  condition  that 
the  linear  response  of  the  particle  number  and  en- 
ergy to  a long- w.ii'elength  disturbanci*  is  zero,  thus 
ensuring  the  appropruite  conservation  laws.  If  we 
use  .in  approximate  G.  then  we  must  construct  a 
A that  maintains  Kq.  (18)  and  this  .illows  us  to  re- 
l.ile  <;  to  the  vertex  lunction  in  the  following  way: 
We  let  ■ — . j.  then  Kq.  (18)  becomes 


1 (’  z ) 


Die  configur.ition  averaged  resolvent  can  be  writ- 
ten 

' I’,- /’-.\,y.',)r' . (201 

lAhu’h  loads  ihe  ('([uation 


i k4  \ X / 


,/z.  .V  2- ,/(,7,„(.:-,)  \ ,1,^  ) vZ.  i’- 


^ k j ^ I ' 

rf> 

In  we  let  i//  - / and  sum  over  all  / with  , - - 


Once  we  recognize 


A,,(z,.z,).  ^ 


we  find  tli.il  the  vertex  must  satisfy  Ihe  equ.ition 

We  note  that  not  onlv  must  this  rel.ilionship  hold 
for  the  exact  verti'x  function  .ind  self-energy  but 
.ilso  111  any  aiiproxiiiiation  in  which  it  is  desired 
lli.it  the  lwo-|).irticle  correl.ilion  function  satisly 
Ihe  W.ird  identilv.  Aq.  (18).  This  mi-.ins  (hat  once 
in  .ipproximalion  is  made  to  \*.  we  c.in  di'lcr- 
miiie  transport  functions  that  .illow  the  consen.i. 
lion  laws  to  be  obeyed.  We  could  imil'i'  farttier 
.qiproxiiiiations  to  {(.  ')  tint  we  would  then  have  no 


giiar.intee  ol  conserving  charge,  energy,  etc.  I'lie 
Ward  identity  is  useful  to  generate  ,i  unique  com- 
(latible  vertex  lunction  only  when  the  frequencies 
are  the  same.  For  tlie  case  ,il  hand.  A;^(C|) 

' a fimction  of  Ihe  sile-di.igon.il  ,iver- 

aged  Green's  function  so  the  vortex  is 

- (»  z ) - (■  is  !,  ~l-  - t c>r,\ 

-ir,  fa'  1 • I ' Y j "tu  , .«  ,)  • (Z.)l 

lot'  I ' . w e make  the  .ip|i ro.v iiiiat  ion 


This  IS  cerl.iinlv  coiisislent  wilti  llii'  W ire!  iih-nliU  . 
•ind  furthermore,  u .illows  us  to  h,  v.  tliat  coniri 
liiitions  to  the  i oiuhel  iv  ity  Iroiii  the  leilex  conic 
lioii.s  llien  mil'll,  the  t w o- p.i  ri  o' K v'trri'l.ilmn 
lunction  i.s  now 
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To  find  (he  I’onductivity  wt*  need  y(6(',  - //)/>, fi(\ , - //)/>_.),  or 

0)/>?,-«f !«(',-//) I,/;). 


(27f 


(2«) 


which  requires 


i;/rrt  Oim 

Tlie  second  term  breaks  up  into 


ilk  jmr 


(30) 


Lei  us  transform  the  W'annier  sum  in  .t  to  a Bloch  sum.  Then  since  (k  |/> , J k')  = /// T (k)fijf, . and  C,  is  diat;- 
onal  in  the  Bloch  representation. 


i HI 


_ 1 
ej  A_- A-(A*)-ej  • 


(31) 


The  propagators  are  even  under  inversion  (k--  k)  Init  the  velocity  T^=  Jt  ik , is  odd,  givint;  us  zero,  ,ind 
all  verte.x  corrections  now  vanish.  In  this  case,  the  fortunate  cancellation  of  vertex  corrections  comes 
about  as  a consequence  of  the  approximation  of  the  proper  self-energy  by  a site- diagonal  quantity,  Kq.  (12), 


V II  /I  KO  1 1 \ll>|  K V 1 1 Rl  iK  ( ()M)l  ( 1 IVI I Y 

We  are  now  in  a position  to  evaluate  u(0).  Because  tlie  vertex  corrections  vanisli, 

/„(' , V)  - Y.  - H)  I 6(A  - //)  j,/»  = HI-  Y l',(k)l’;(iT){<k|  6(A  - //) , k))  '.  (32) 

ijlm  f 


With  the  definition 

!;k,fi{A-//)  k)>  . (-1  i7)lm(7jf(A-), 
we  get  for  f - 2 |<r(j.')  = OiJo.')  | , 


X Y '■.(k)"|Ini(T;i:('‘)r.  (33) 

wticre  we  have  included  a factor  of  2 for  the  two 
IHissilile  spin  orientations.  ,-\t  / 0,  ■' 

'('  - . ),  with  .1  the  chemical  potential,  and  the 
conduettvits  [n'r  atom  is 

r(0)  Y l',lk)’|lm(Ti;f(M')r  (;)4) 

or 


( • (/i*)  - A |‘  f I Ml  A * (^E)  y 


1 ■ 15.  r!i>  !i  < !»  « • bi.lii.  ii  1 it \ 

-Itnu-MMonril  .-ilInN  u Ith  I j .p,l  ,•  ji.j.  u-n.*  It\ 

■ij.'Uc  > irom  1 i.-..  V hnwn  in  th<-  n.iM.I  rm  i 

:;\  l^  in  unil>  >>1  h:ilf  5trin'’.v  liili. 
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This  natural  separation,  only  ixjssible  tor  a proiier 
self-energy  indepeiulent  of  isolates  the  lifetime 
and  energy  shifts  of  the  sinule-partiele  exeilations 
from  that  part  of  the  conductivity  which  pei  tains 
to  the  particular  lattice  under  study.  We  will  con- 
centrate on  a three-dimensional  simple  cubic  lat- 
tice wiUi 

',(cos/.\+  cosA’,+  cos/,’^1 . 

Then, 

r,(k)-  = ^sin'/.-^  = .L(i  _ co.s'V.-^. 


Consider  the  functions 

xl-riT:- 


for  which  w'e  have 


V. 

The  imaiiinary  part  of  both  P.IK')  and 
vanish  outside  the  unperturbed  band  and 


■J(01 


1mA’  ( 4*) 


( 4 - [?e.\’  (4*)  - E I"  + lm.\ 


(30) 


We  have  calculated  the  dc  conductivity  for  our 
three-dimensional  alloy  in  order  to  illustrate  our 
formal  results,  i^lenerally.  there  are  two  ways  in 
which  the  dc  conductivity  can  vanish.  If  the  den- 
sity of  states  at  the  Fermi  level  is  zero  then  so  is 
o(0).  In  addition,  we  can  have  a finite  p{4),  but  a 
zero  mobility  because  of  wave-function  localization. 
Equation  (36)  only  admits  a zero  in  i(0)  if  /dp)  is 
zero  so  we  cannot  t.ike  the  latter  possibility  into 


account.  The  conductivity  is  displ.ived  m Fit;,  II 
at;ain.st  its  respective  density  of  st.ites.  I'here  is 
a rather  direct  correlation  between  the  niapnitude 
of  the  density  of  states  and  that  of  the  conductsviiv 
This  relationship  is  understood  in  terms  ol  the 
avail.iliility  of  states  at  a itiven  eneruv  to  which  .in 
initial  .slate  can  make  a Ir.insition  Ue  also  find 
peaks  in  r(0)  which  we  associ.ite  with  velociu 
peaks  in  the  cubic  band  structure. 
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